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ABSTRACT 

We model the infrared (IR) emission from the double-ring disk of HD 141569A, 
using a porous dust model that was previously shown successful in reproducing the 
spectral energy distributions (SEDs) of the disks around (3 Pictoris and HR4796A. 
The dust consists of either unaltered or highly processed interstellar materials and 
vacuum with a volume fraction of ~ 90%. Together with a population of polycyclic 
aromatic hydrocarbon (PAH) molecules of which a large fraction is charged, the porous 
dust model made of coagulated but otherwise unaltered protostellar interstellar grains 
provides an excellent fit to the entire SED from the mid-IR to millimeter wavelengths 
including the PAH mid-IR emission features. The nondetection of the 21 ^m crystalline 
silicate feature predicted from the porous dust model composed of highly-processed 
grains (mainly crystalline silicate dust) imposes an upper limit of ~ 10% on the mass 
fraction of crystalline silicates. Spectroscopic and broadband photometric predictions 
are made for SIRTF observations which will provide further tests of the applicability of 
the porous dust model to circumstellar disks. 

Subject headings: circumstellar matter — dust, extinction — infrared: stars — plane- 
tary systems: protoplanetary disks — stars: individual (HD 141569A) 



1. Introduction 

Since dust disks around young stars, depending on their age, are the source material or the 
remnants of newly- formed planets, the physical, chemical, and dynamical properties of circumstellar 
disks and their constituent grains are crucial in understanding the formation process of planetary 
systems. 

Over the past 2 decades, impressive evidence has been assembled for the commonality of 
circumstellar dust disks around pre- main-sequence stars (T Tauri stars and Herbig Ae/Be stars), 
main sequence (MS) stars, post-MS stars (red giants), and a white dwarf (see Backman &: Paresce 
1993, Habing et al. 2001, Zuckerman 2001, Lagrange, Backman, & Artymowicz 2001 for reviews). 
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To date, at least 15% of the A-K type stars are found to have dust disks around them (see Lagrange, 
Backman, & Artymowicz 2001). The majority of these disks are spatially extensive, with their dust 
quantities and intensities of scattered-starlight and of dust thermal emission gradually decreasing 
outward from the central stars. However, so far there exist four exceptions: (1) e Eridani - a 
~ 800 Myr old (Henry et al. 1996) K2V MS star - has a nearly face-on dust ring with a mass and 
radius (r ~ 60 AU) similar to the solar system Kuiper Belt and is thought to be an analog to the 
young solar system (Greaves et al. 1998); (2) Formahault (a PsA) - a ~ 200 Myr old (Barrado y 
Navascues et al. 1997) A3V MS star - exhibits a "doughnut" -like dust annulus with a radial distance 
of ~ 150 AU from the star and a width of ~ 60 AU (Holland et al. 1998, 2003); (3) HR4796A - 
a ~ 8 ± 3 Myr old A0V star (see Jura et al. 1998) - displays a ring-like disk peaking at ~ 70 AU 
from the central star and abruptly truncated both inward and outward with a width of < 17 AU 
(Schneider et al. 1999; Koerner et al. 1998; Jayawardhana et al. 1998; Telesco et al. 2000); (4) 
HD 141569A - a nearby (distance to the Earth d ~ 99 pc; van den Ancker, de Winter, & Tjin A 
Djie 1998) ~ 5 ±3 Myr old (Weinberger et al. 2000) pre-MS Herbig Ae/Be star with a spectral type 
of B9.5V, an effective temperature of T e g- « 10000 K and a surface gravity of Igg « 4.0 (Dunkin, 
Barlow, k Ryan 1997). 

The presence of dust around HD 141569A was first reported by Jaschek, Jaschek, & Egret 
(1986), de Grijp, Miley, & Lub (1987), and Walker & Wolstencroft (1988) based on the IRAS {In- 
frared Astronomical Satellite) detection of infrared (IR) excess over expected from its photosphere. 
It has recently been imaged in scattered light and in thermal IR emission. With the discovery of its 
large and morphologically complicated debris disk (see below), HD141569A has recently aroused 
considerable interest. 

The disk around HD 141569A appears to have two ring-like structures: an outer ring at ~ 
325 AU from the star with a width of ~ 150 AU; and an inner ring at ~ 185 AU. The outer ring 
was reported independently by Augereau et al. (1999a) and by Weinberger et al. (1999), respectively 
based on the 1.6 /im and 1.1 /im imaging of the scattered- light. They both used the Near-IR Camera 
and Multiobject Spectrometer (NICMOS) on the Hubble Space Telescope (HST). The inner ring was 
discovered by Weinberger et al. (1999) at 1.1 /im using the same instrument. Weinberger et al. 
(1999) also reported the detection of a gap or "dip" (a region of depleted material) at ~ 250 AU 
from the star with a width of ~ 60 AU. They interpreted the gap as a sign of one or more planets 
orbiting within the disk (the planets do not have to be in the gap), but other explanations such as 
dust migration also exist (e.g. see Takeuchi & Artymowicz 2001). 

This ring-gap-ring structure was later confirmed by optical coronagraphic observations made 
by Mouillet et al. (2001) at an effective wavelength of 5850 A using the HST Space Telescope Imaging 
Spectrograph (STIS). With a higher spatial resolution and a higher signal-to-noise ratio than the 
NICMOS/HST data, Mouillet et al. (2001) found that the inner ring, with a FWHM (full width 
of half maximum) of ~ 50 AU, actually peaks at ~ 200 ± 5 AU. Strong brightness asymmetry was 
also reported by Mouillet et al. (2001). This asymmetry, interpreted as a result of the gravitational 
perturbation caused by putative massive planets within the disk, by the two stellar companions 
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HD141569B and HD141569C, or both types of objects (Mouillet et al. 2001), was confirmed by 
the ground-based 2.2 /im near-IR imaging observation at the Palomar 200 inch telescope using the 
adaptive optics system (Boccaletti et al. 2003). 

In addition to the two annuli, the existence of a third dust component - a population of warm, 
mid-IR emitting dust within ~ 100 AU from the star - was recently reported by Fisher et al. 
(2000). It was shown that the 10.8 and 18.2 /im emission observed with the University of Florida 
Observatory Spectrometer and Camera for the Infrared (OSCIR) on the Keck II telescope extend 
out to ~ 100 AU from the star (Fisher et al. 2000). The existence of such an extended, warm source 
close to the star, originally suggested by Augereau et al. (1998) to account for the disk IR emission 
at wavelengths A < 10 fim, was later confirmed by the 12.5, 17.9, and 20.8 /im imaging using the 
JPL Mid-InfraRed Large-well Imager (MIRLIN) on the Keck II telescope (Marsh et al. 2002). 

More recently, the HRC (High Resolution Channel) coronagraphic B, V, and I images of the 
HD 141569A disk obtained with the HST Advanced Camera for Surveys (ACS) revealed a tightly- 
wound spiral structure with two arms, one of which stretches outward to reach the nearby com- 
panion stars (Clampin et al. 2003). 

In a previous paper (Li & Lunine 2003), we have shown that the IR emission from the HR4796A 
disk can be well fit by a model invoking highly porous cometary-type dust. The porous cometary 
dust model was originally proposed for the (3 Pictoris disk and was shown successful in reproducing 
its spectral energy distribution from the near-IR to millimeter wavelengths including the 10 /im 
amorphous and the 11.3 /im crystalline silicate features (Li & Greenberg 1998). In this paper we 
will model the observed IR and millimeter photometric and mid-IR spectroscopic signatures of 
HD 141569A. The objectives of this paper are three-fold: (1) We wish to infer the physical and 
chemical properties of the dust in the HR4796A disk and its relationship to the formation of planets. 
(2) We wish to know how widely applicable is the porous dust model; if it can be shown that the 
porous dust model is also valid for many other disks at different evolutionary stages and with 
different geometrical structures, it will be a valuable guide for interpreting future data sets. This 
is of particular importance given that the incoming Space Infrared Telescope Facilities (SIRTF) 
will accumulate a rich set of IR photometric and spectroscopic data for a large number of dust 
disks; we would hope that the porous dust model presented in this and related papers could serve 
as a starting point in analyzing these data. (3) We wish to make near- and mid-IR spectral and 
broadband photometry predictions for the HD 141569A disk; these predictions can be compared to 
future SIRTF observations in order to test further the validity of the porous dust model. 

This paper is organized as follows: we first summarize in §2 the photometric and spectroscopic 
data available for the HD 141569A disk that are relevant in modeling the IR emission of the disk. 
We then give an overview discussion in §3.1 of the proposed dust model with emphasis on dust 
composition and morphology from the point of view of protoplanetary disk evolution, followed by 
detailed discussions of the PAH component in §3.1.1 and the porous dust component in §3.1.2. The 
IR emission modeling technique is presented in §3.1.3. In §3.1.4 the PAH photoelectric emission 
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rates and the electron recombination rates are calculated, using the methods described in §A and 
§B. The key model parameters are discussed and specified in §3.1.5. Model results are presented 
in §3.2. In §4 we discuss (1) previous studies of the SED modeling of the HD141569A disk, 
(2) comparison of the HD 141569A disk with the HR4796A disk, (3) exploration of parameter 
space (§4.1), (4) PAH photodestruction (§4.2), and (5) dust loss through radiation expulsion and 
Poynting- Robertson drag (§4.3). In §4.4 we calculate the dust IR flux densities integrated over the 
SIRTF/IRAC and MIPS bands predicted for our best-fitting models. Spectroscopic predictions are 
also made for SIRTF/IRS observations. The major conclusions are summarized in §5. 



2. Available Photometric and Spectroscopic Data 

To facilitate modeling of the dust properties of the HD 141569A disk, we have assembled the 
available optical, IR, and millimeter photometric data and mid-IR spectra. These include (1) the 
BV optical photometry (Sylvester et al. 1996); (2) the JHKLL'M near-IR photometry obtained by 
Sylvester et al. (1996) using the single-channel bolometer UKT9 at the United Kingdom Infrared 
Telescope (UKIRT); (3) the 12, 25, 60, and 100 fim IRAS photometry (Sylvester et al. 1996); 
(4) the 10.8 and 18.2 fim OSCIR/Keckll imaging by Fisher et al. (2000); (5) the 12.5, 17.9, and 
20.8 fim MIRLIN/Keckll imaging by Marsh et al. (2002); (6) the 7.5-13.5//m (spectral resolution 
AA = 0.17/Ltm) and 15.8-23.9 fim (AA = 0.27/mi) spectroscopy obtained by Sylvester et al. (1996) 
using the UKIRT Cooled Grating Spectrometer 3 (CGS3); (7) the 1350 [xm SCUBA (Submillimetre 
Common-User Bolometer Array) measurement by Sylvester, Dunkin, & Barlow (2001) using the 
James Clerk Maxwell Telescope (JCMT). Taking E{B - V) = 0.095 (Weinberger et al. 1999), 
we have corrected the interstellar reddening assuming a Ry = 3.1 (total-to-selective extinction) 
interstellar extinction law (see Table 6 in Li & Draine 2001). 

Table 1 summarizes the photometric data listed above, the instruments/telescopes and their 
beam sizes through which the observational data were obtained, and the stellar atmospheric and 
dust disk contributions to each photometric band. In Figure la we plot the reddening-corrected 
optical, IR, and millimeter photometric results as well as the T e g = 10000 K and Igg = 4.0 Kurucz 
model atmospheric spectrum. Figure lb illustrates the CGS3 mid-IR spectrum which clearly ex- 
hibits the 11.3 /xm out-of-plane C-H bending feature and the long- wavelength wing of the 7.7//m 
C-C stretching feature (see Li Sz Draine 2001a for the PAH band assignment). For the sake of 
comparison, the 10.8, 18.2 /im OSCIR/Keckll data, the 12, 25 p IRAS data, and the 12.5, 17.9, 
20.8 fim MIRLIN/Keckll data, are also plotted in Figure lb. It is seen that these mid-IR data 
closely agree with one another. 1 The largest discrepancy (f» 30%) occurs for the 17.9 /jm MIR- 
LIN/Keckll (Marsh et al. 2002) data and the 18.2 /im OSCIR/Keckll data (Fisher et al. 2000). 
This discrepancy cannot be explained in terms of the difference in the aperture sizes of these two 



lr The 12 fim IRAS data is about 12% higher than the band-averaged CGS3 spectrum. We have therefore reduced 
the former by a factor of 1.12 to bring the 12 pm IRAS data into agreement with the CGS3 spectrum. 
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instruments. As a matter of fact, the former (which has a higher 17.9 //m flux density) was obtained 
using a smaller aperture. The CGS3 mid-IR spectrum, obtained with a 5.5" beam, appears to be 
intermediate between the MIRLIN and the OSCIR data. Therefore, we ascribe this discrepancy to 
observational uncertainties. Both data points will be included in our detailed dust modeling. 



Table 1: Existing optical, IR and millimeter photometric data for the HD 141569A system (see §2). a 



A 


F obsb 


jp st cir c 


F dustd 


F obab 




F duatd 


Telescope/ 


Beam/ Aperture 


(Mm) 


(Jy) e 


(Jy) 


(Jy) 


(erg s 1 cm 2 (im 1 ) 


(erg s 1 cm 2 (im 1 ) 


(ergs" 1 cm - 2 tim -1 ) 


Instrument 


Size 


0.45 (B) 


10.50 


10.50 




1.56 X 10~' 


1.56 X 10~ ' 








0.55 (V) 


6.70 


6.70 




6.64 X 10~ S 


6.64 x 10~ s 








1.22 (J) 


3.19 


3.19 




6.43 X 10~ 9 


6.43 X 10 -9 




UKIRT 


7.8" 


1.65 (H) 


2.02 


2.02 




2.23 X 10~ 9 


2.23 X 10~ 9 




UKIRT 


7.8" 


2.18 (K) 


1.29 


1.29 




8.14 x 10~ 10 


8.14 x 10~ 10 




UKIRT 


7.8" 


3.55 (LJ 
3.76 (I/) 


0.59 


0.58 


0.018 


1.41 x 10~ 10 


1.37 X 10~ 10 


4.29 x 10~ 12 


UKIRT 


7.8" 


0.54 


0.50 


0.040 


1.14 x 10~ 10 


1.06 X 10~ 10 


8.54 x 10 -12 


UKIRT 


7.8" 


4.77 (M) 


0.40 


0.33 


0.070 


5.33 X 10 -11 


4.41 x 10 -11 


9.25 x 10 -12 


UKIRT 


5.0" 


12 


0.49 


0.058 


0.43 


1.02 x 10" 11 


1.21 x 10~ 12 


9.02 x 10 -12 


IRAS 


1' 


25 


1.87 


0.014 


1.86 


8.98 x 10~ 12 


6.64 x 10~ 14 


8.91 x 10~ 12 


IRAS 


1' 


60 


5.54 


0.0024 


5.54 


4.62 x 10~ 12 


2.04 x 10~ 15 


4.61 X 10~ 12 


IRAS 


1' 


100 


3.48 


0.0009 


3.48 


1.04 x 10~ 12 


2.65 X 10~ 16 


1.04 x 10^ 12 


IRAS 


1' 


10.8 


0.318±0.016 


0.071 


0.25 


8.18 X 10~ 12 


1.83 X 10~ 12 


6.34 x 10~ 12 


OSCIR/Kcckll 


7.9" x 7.9" 


18.2 


0.646±0.035 


0.026 


0.62 


5.85 X 10 -12 


2.34 x 10 -13 


5.62 x 10" 12 


OSCIR/Kcckll 


7.9" x 7.9" 


12.5* 


0.333±0.022 


0.054 


0.28 


6.39 X 10~ 12 


1.03 x 10~ 12 


5.36 X 10" 12 


MIRLIN/Keckll 


5" 


17. 9 f 


0.936±0.094 


0.027 


0.91 


8.76 X 10 -12 


2.50 x 10~ 13 


8.51 X 10 -12 


MIRLIN/Kcckll 


5" 


20. 8 f 


1.19±0.16 


0.020 


1.17 


8.25 X 10~ 12 


1.38 x 10~ 13 


8.11 X 10~ 12 


MIRLIN/Keckll 


5" 


1350 


0.0054 ± 0.001 





0.0054 


8.89 X 10 _ls 


8.03 x 10 -21 


8.88 X U)- lx 


SCUBA 


22" 



a In the SED modeling presented in later sections, we will include both the CGS3 mid-IR spectrum and the IRAS, 
OSCIR/Keckll, MIRLIN/Keckll, and SCUBA photometry. 
6 Total flux density observed from the Earth. 
c Stellar photospheric contribution to the observed flux density. 
d Dust IR excess from the disk. 




The MIRLIN/Keckll instrument has a field of view of 18" x 18". But the flux given in this table was obtained by an 
integration over a 5" diameter circular aperture (see Marsh et al. 2002). 



3. Modeling the Dust IR Emission 
3.1. Dust Model 

Interstellar grains undergo coagulational growth in cold, dense molecular clouds and protostel- 
lar and protoplanetary dust disks. This process creates dust with fluffy, inhomogeneous structures 
and fairly quickly leads to the formation of planetesimals of cometary to asteroidal sizes. 

Based on its short lifetime due to radiative expulsion, the dust in the HD 141569A disk is 
thought to be of secondary debris origin (Weinberger et al. 1999, 2000; Fisher et al. 2000; Marsh 
et al. 2002); i.e, the dust is generated by the breaking down of large bodies like planetesimals, 
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comets, and asteroids, rather than direct remnant left over from the star formation process. Since 
planetesimals and comets are formed through coagulation of porous dust aggregates (Cameron 
1975, 1995; Weidenschilling & Cuzzi 1993; Greenberg & Li 1999; Beckwith, Henning, & Nakagawa 
2000), the dust generated by collisions of planetesimals should resemble the original dust from 
which they are built up. 

Following Li & Lunine (2003), we will model the HD 141569A dust as porous aggregates of 
either processed or unprocessed interstellar materials. We will consider two dust types - each 
represents an extreme case - the "cold-coagulation" dust model and the "hot-nebula" dust model. 
While the former assumes that the dust aggregate is formed through cold-coagulation of unaltered 
protostellar interstellar grains, the latter goes to the other extreme by assuming that the constituent 
grains in the dust aggregate are highly-processed in the protostellar nebula where the originally 
amorphous silicate grains all are annealed and the carbon grains all are destroyed by oxidization 
and are converted into CO (see Li & Lunine 2003 for details). 

A porous dust grain is characterized by (1) its fluffiness or porosity P (the fractional volume 
of vacuum); (2) its spherical radius a (the radius of the sphere encompassing the entire aggregate; 
we assume all grains are spherical in shape); (3) composition - the "cold-coagulation" dust model 
assumes the dust to be composed of amorphous silicate and carbonaceous materials (and H2O- 
dominated ices in regions colder than ~ 110 — 120 K); the "hot-nebula" dust model assumes the 
dust to be composed of only crystalline silicates (and ices in cold regions). As discussed in Li 
& Lunine (2003), for the "cold-coagulation" dust model, the mixing mass ratios for the silicate, 
carbon and ice constituent grains are approximately derived, by assuming cosmic abundance, to 
be mcarb/w-sii w 0.7 and mi cc /(m s ii + m cav ^) « 0.8 where m s n, m car b, and m- iCC are respectively the 
total mass of the silicate, carbon, and ice subgrains. For the "hot-nebula" dust model, we will take 
m ca,rb / m s\\ = and m\ ce /m s \\ = 1/3 (see Li & Lunine 2003). We adopt a mass density of p s ji = 
3.5 gem -3 , p car b = 1.8 gem" 3 , and p\ ce = 1.2 gem -3 respectively for the silicate, carbonaceous, and 
ice materials. 2 

In addition to the porous dust component discussed above, the presence of PAH molecules in 
the HD141569A disk is clearly indicated by the red wing of the 7.7 //m feature and the 11.3 fim 
feature seen in the CGS3 mid-IR spectrum (Sylvester et al. 1996; see Fig. lb). Therefore, we will 
consider two dust populations in our models: the porous dust component (§3.1.1) and the PAH 
component (§3.1.2). 



2 For pure H2O ice we have pi cc w l.Ogcm -3 . Since interstellar ice is a "dirty" mixture of H2O, NH3, CO, CO2, 
CH3OH, CH4 and many other molecules (see Appendix A in Li & Lunine 2003), we adopt pi ce = 1.2 gem -3 for "dirty 
ice" material (Greenberg 1968). 
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3.1.1. Poly cyclic Aromatic Hydrocarbon Dust 

PAHs are excited by ultraviolet (UV), visible, and to a lesser extent by long wavelength photons 
(Li & Draine 2002a). Following Li & Draine (2002a), we define U(r) as the 912 A - 1 fim starlight 
intensity incident on the HD 141569A disk at a distance of r from the central star relative to the 
value for the Mathis, Mezger, & Panagia (1983; hereafter MMP) solar neighborhood interstellar 
radiation field 3 

(«./»)• J^-^ ./ 4.75x 10* AU \' 

J 912A A 

where 7.32 x 10" 3 AU) is the HD 141569A stellar radius; F£ is the flux per unit wavelength 

(ergs -1 cm -2 /iin -1 ) at the top of the illuminating star's atmosphere which is approximated by the 
Kurucz model atmospheric spectrum for B9.5 V stars (T e g = 10000 K and Igg = 4.0; Kurucz 1979); 
c is the speed of light; uf RF is the energy density of the MMP interstellar radiation field. In Figure 
lc we plot the starlight spectrum at r = 4.75 x 10 4 AU which corresponds to U = 1. Since the 
HD 141569A disk extends to a distance of only ~ 500 AU, the dust in this disk receives far more 
intense irradiation than the dust in the diffuse interstellar medium: U > 9000! 

Illuminated by starlight with an intensity of U, the IR emissivity per gram dust (erg s _1 /im" 1 g -1 } 
for a mixture of neutral and ionized PAHs with a size distribution of dnpAn/da is 

£ AH (U) = da^ f°° dT 4itB\(T) 

J a PAH da J 

min w 

x{0 ion (^,a)C a ^ H+ (a, A) ( \ + [1 - ion (C7, a)]C a p ^ H (a, A) f \ } 

dnpAH47T 3 l-i 

da— r-a°p PA H \ (2) 

pah da 6 > 

min 

where a is the spherical radius of a PAH molecule; 4 is the lower-cutoff PAH size; C^ u ° 

and C a ^ H+ are the absorption cross sections for neutral and ionized PAH molecules, respectively 
(see Li & Draine 2001a); B X (T) is the Planck function at temperature T; dP PAn ° (U,a,T) and 
dP PAR+ (U,a,T), are, respectively, the probabilities that the vibrational temperature will be in 
[T,T + dT] for neutral and charged PAHs illuminated by starlight intensity U; (pi on (U,a) is the 
probability of finding a PAH molecule of radius a in a non-zero charge state; />pah is the mass 
density of PAHs which is taken to be that of graphite (« 2.24 gcm~ 3 ). 



f 



3 For radiation fields with the same 912 A - 1 /im intensity, the 6-13.6 eV starlight intensity of the T e g = 10000 K 
Kurucz radiation field is about 1.27 times that of the MMP field. 

4 Small PAHs (with < 100 carbon atoms) are expected to be planar (see Appendix A in Draine & Li 2001). The 
term "PAH radius" used in this paper refers to the radius a of a spherical grain with the same carbon density as 
graphite (2.24 gem -3 ) and containing the same number of carbon atoms A^c: a = 1.288N^/ 3 A [or At = 0.468(a/A) 3 ]. 
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Based on the "thermal-discrete" method (Draine & Li 2001), the temperature distribution 
functions dP/dlnT have been computed by Li & Draine (2002a) for neutral and charged PAHs 
and ultrasmall silicate grains of a range of sizes 3.5 A < a < 350 A illuminated by stars with a 
range of effective temperatures 3000 K < T e ff < 30000 K and with a range of starlight intensities 
0-1 < U < 10 5 . In this work we take the Li & Draine (2002a) dP/d\nT results for T cff = 10000 K. 

We adopt a log-normal size distribution for the PAHs in the HD 141569A disk, characterized 
by two parameters: ao and a; ao and a respectively determine the peak location and the width of 
the log-normal distribution: 



dnpAu/da 




ln(a/a ) 12 



^72a{l-erf[ln(aPAH /ao)/V 2a]}a-^ ~ " for ° > ^ 

(3) 

where erf(y) = J Q y exp (— x 2 ) dx is the error function. The log-normal functional form for the PAH 
size distribution was shown successful in modeling the PAH mid-IR emission spectra observed for 
the Milky Way diffuse interstellar medium (Li & Draine 2001a), reflection nebulae (Li <fe Draine 
2002a), and the Small Magellanic Cloud (Li & Draine 2002b). 

Lacking a priori knowledge of the balance between the photodestruction and collisional pro- 
duction of PAHs as a function of PAH size in the HD 141569A disk, we adopt a lower-cutoff of 
a mhi 1 = 3.5 A (corresponding to Nq ~ 20 for PAHs) which is the minimum survival size for PAHs 
in the diffuse interstellar medium (see Li & Draine 2001a). It will be seen in §4.2 that PAHs smaller 
than ~ 4.6 A will be photolytically unstable in the inner r < 100 AU region around the star during 
the lifetime of the disk. PAHs are also subject to radiative expulsion. Continuous replenishment 
of this material, presumably by sublimation of the icy mantles (in which interstellar PAHs have 
condensed during the dense molecular cloud phase) coated on the individual subgrains of large 
porous dust produced by collisions of large parent bodies, is required to maintain the PAH disk. 
Therefore, it is reasonable to assume that small PAH molecules with a < 4.6 A are also continu- 
ously replenished so that there exists a stable distribution of small PAHs in the disk. Since the dust 
(including PAHs) in the disk originates from the interstellar medium (ISM), it is thus reasonable 
to adopt the interstellar lower-cutoff aj^ 1 = 3.5 A for the PAHs in the disk. But we will also see 
in §4.2 that models with = 4.6 A are also capable of closely reproducing the observed SED 

including the PAH emission features. 



3.1.2. Porous Dust 

Large porous dust in the HD 141569A disk will attain a steady-state temperature determined 
by balancing absorption and emission, 

^ CT s (a,X)F^dX = J o CS°(a,\)4xB x (T[a,r])d\ (4) 
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where C^(a, A) is the absorption cross section of porous dust of spherical radius a at wavelength 
A; T(a, r) is the equilibrium temperature of dust of size a at a radial distance of r from the star. 

The IR emissivity (ergs -1 nmT 1 ) for porous dust with a size distribution of dn-p^/da located 
at a distance of r from the central star is 



JA PD (r) 



/ da -^4nB x (T[a,r})C™(a,\) . (5) 



We assume a power-law dust size distribution for the porous dust 



da In [a max /a 

(l-a 



a , a = 1 

minj 



l-a _ l-a 
"max « min 



a ^ 1 ; (6) 



where a is the spherical radius, a m i n is the lower-cutoff, a max is the upper-cutoff, and a is the 
power-law index. We take a m i n = 1 |im, 5 and a max = 1 cm (this is not a critical parameter since 
grains larger than ~ 100 fim emit like blackbodies and their IR emission spectra are size-insensitive; 
see Figure 6). 

We use Mie theory to calculate the absorption cross sections for porous dust. The fluffy 
heterogeneous dust aggregate is represented by an equivalent homogeneous sphere with an effective 
dielectric function e e ff • Dielectric functions for the constituent dust materials are taken from (1) 
Draine & Lee (1984) for amorphous silicate dust; (2) Li & Draine (2001b) for crystalline silicate dust; 
(3) Li & Greenberg (1997) for carbonaceous dust; (4) Li & Greenberg (1998) for H^O-dominated 
ice. Let e s n, e car b, and £i ce be the complex dielectric functions of silicate (either amorphous or 
crystalline), carbonaceous, and ice dust, respectively. Let / s n, / C arb ; and /ice respectively be the 
volume fraction of silicate, carbonaceous, and ice dust in an aggregate, which can be obtained from 
the mass mixing ratio m carb /m si i and m icc / (m sU + m carb ) and the fluffiness P (= l-/ s ii-/ ca rb-/ice) 
of the aggregate (see §3.1). We first employ the Maxwell-Garnett effective medium theory (Bohren 
& Huffman 1983) to calculate the average dielectric functions for the ice-coated silicate subgrains 
e' sil and the ice-coated carbonaceous subgrains e' carh 

I _ (! + 2 /sil) £ sil + 2(1- /sil) g icc f>— fl(fA-Cf\. <7\ 

fc sil — e ice / „, \ | /(-j . p/ \ i J sil — Jsil/ Usil 1" CJiccJ , {< ) 

(1 - / S ii) esii + (2 + / s ' n ) e ice 

,/ _ , (j; j~ ^^carb) g carb + 2 (1 — / carb ) 6jcc w _ , , r , , /irtf 1 , (o\ 

fc carb — fc ice / \ /„ \ > Jcarb — Jcarb/ [Jcarb T \± C,) Jice\ , \0) 

I 1 _ /carb/ £ carb + + J C arbJ £ ice 

where £ is the fraction of ice condensed on the silicate cores. For the "cold-coagulation" dust model, 
we assume £ = 0.5 (i.e. in an dust aggregate, the available ice equally condenses on the silicate 



5 A fluffy grain of a = 1 /im with a porosity of P = 0.90 consists of ~ 100 constituent individual (interstellar) 
particles which have a typical size of a ~ 0.1 (im (see Li & Greenberg 1997). Models with a smaller a m i n (= 0.1 /im) 
and a larger a m i n (= 10 /jm) will be discussed in §4.1. 
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subgrains and on the carbonaceous subgrains); for the "hot-nebula" dust model, we assume £ = 1 
since we assume that all carbon dust has been destroyed by oxidation (see Li & Lunine 2003). We 
then use the Bruggman effective medium theory (Bohren & Huffman 1983) to calculate the mean 
dielectric functions e e s for the porous heterogeneous dust aggregate 

(/sil + £/icc) ~T~ I o l~ [/carb + (1 — £) /ice] ~j ; X ^ ~~ /sil ~~ /carb ~~ /ice) ; ^ = ; 
£ sil + 2£ cff e carb + 2e eS 1 + 2e cS 

(9) 

where the third term accounts for the contribution from the vacuum (e = 1). By definition, Eqs.(7- 
9) are valid for both the "cold-coagulation" dust model (£ = 0.5) and the "hot-nebula" dust model 
(£ = 1) either with (/; ce > 0) or without (/; ce = 0) ice mantles. 



3.1.3. IR Emission Spectrum 

For a given dust size distribution and a given disk structure (i.e., the dust spatial density 
distribution), the emergent IR emission spectrum can be obtained by integrating over the dust size 
range, and over the entire disk. Assuming an optically thin (see §3.2) and radially symmetric disk, 
the flux density (ergs -1 cm -2 /im" 1 ) received at the Earth is 



1 



4vrd 2 



jT{r)o{r) 27rrdr + jl A}1 (U)m PA}i (U) 



(10) 



where d ~ 99 pc is the distance from the star to the Earth; a(r) is the surface density for the 
porous dust component; r- m and r out are respectively the inner and outer boundaries of the disk; 
w-pah(^) is the total mass of the PAH dust required to account for the observed PAH IR emission 
bands if the PAH molecules in the HD 141569A disk are illuminated by starlight of an intensity of 
U times that of the 912 A - 1 /im MMP interstellar radiation field. We will justify the optical thin 
treatment in §3.2 (also see Table 2). The HD 141569A disk is actually asymmetric (Mouillet et al. 
2001; Boccaletti et al. 2003; Clampin et al. 2003), therefore, our model, based on an assumption of 
a radially symmetric disk, is somewhat simplified. But this is not expected to affect the conclusion 
of this paper since the asymmetric dust spatial distribution which would account for the observed 
brightness asymmetry, with some fine tuning, would also be able to reproduce the observed SED 
as can be seen from Eq.(10). 

Due to their single-photon heating nature, the IR emission spectral shape of PAHs is indepen- 
dent of the starlight intensity U ; the absolute emissivity level simply scales with U (Draine & Li 
2001; see Fig. 13 in Li & Draine 2001a). Namely, j£ AH {U)/U remains identical over the entire dust 
disk, provided that the charging condition U ^/T gas /n e for PAHs does not vary in the disk. 6 For 
computing the PAH component, it is therefore not necessary to integrate over the disk; instead, it 



6 The PAH ionization fraction <^i on depends on U /n e and T gas where n e is the electron density in the disk and T gas 
is the gas temperature (Weingartner & Draine 2001). 
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is sufficient just to calculate the PAH IR emission spectrum for a single radiation intensity U (i.e., 
a single region r in the disk). The required PAH mass is inversely proportional to U. 

The total mass for the porous dust component is 



f° Ut ( \o A T maX ^ ^P£ 47r 3 nn 
/ a[r) 2ir rar / da— a ppp, (11) 

J r- Jn ■ da o 

" 1 in "mm 



where ppp, is the mass density of the porous dust component (see Appendix B in Li & Lunine 
2003). 

The vertical optical depth is 

TA(r)= / da-^a(r)C™(a,\) (12) 

" ^min 

where C^(a, A) is the extinction cross section of the porous dust of spherical radius a at wavelength 
A. The optical depth from the PAH component is negligible (see §3.2). 

The dust spatial distribution is well constrained by the near-IR imaging of scattered starlight 
(Weinberger et al. 1999; Augereau et al. 1999a) and mid-IR imaging of dust thermal emission (Marsh 
et al. 2002). We adopt the following analytical formula to approximate the dust distribution which 
was derived from the scattered light (see Fig. 3 of Weinberger et al. 1999) and dust thermal emission 
(see Fig. 4 of Marsh et al. 2002) 7 

a(r)/a p = 0.223exp |-41n2 [(r/AU - 55) /70] 2 } , r in < r < 105 AU ; (13) 
= exp{-41n2[(r/AU- 203)/90] 2 } , 105 AU < r < 235 AU ; (14) 
= 0.914 140 2 / (r/AU - 310) 2 + 140 2 ] , 235 AU < r < r out ; (15) 

where <t p is the mid-plane surface density at r = r p = 203 AU. We note that, since we use 
this dust spatial distribution, it is expected that our models are also able to reproduce the imaging 
observations. Eq.(13) represents the 12.5, 17.9, 20.8 /im mid-IR emitting dust component (Marsh et 
al. 2002); Eq.(14) and Eq.(15) describe the two annular structures peaking at r rs 200AU,325AU, 
respectively (Augereau et al. 1999a; Weinberger et al. 1999; Mouillet et al. 2001). The inner 
boundary r; n is set at where grains are heated to > 1500 K. Therefore, r; n is a function of dust 
size. For micron-sized grains, the inner boundary is roughly at r; n = 0.2 AU. The outer boundary 
is taken to be r out = 500 AU which is expected from the disk truncation caused by the tidal effects 
of HD 141569B, a companion star of HD 141569A (Weinberger et al. 2000). 8 This is also consistent 



7 We assume that the dust vertical distribution is determined by the vertical hydrostatic equilibrium (see Eqs.[B2- 
B4] in §B). 

8 It was shown by Artymowicz & Lubow (1994) that circumstellar disks will be truncated by the tidal effects 
of a companion star in circular orbit at approximately 0.9 of the average Roche lobe radius. For the HD 141569A 
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with the near-IR scattered light observations reported by Weinberger et al. (1999) and Augereau 
et al. (1999a) which show that the HD 141569A disk extends to a radial distance of about 500 AU. 

The mean starlight intensity (U), averaged over the different portions of the disk, and weighted 
by the surface density of the porous dust population, is 

f out U(r)a(r)27rrdr 

(U) = 

K 1 [ r ° ut a{r)2Trrdr 

= 1.24 xlO 6 , r in = 1AU < r < 105AU ; 

= 2.76 xlO 4 , 105 < r < r out = 500 AU ; 

= 4.03 x 10 4 , r in = 1 AU < r < r out = 500 AU . (16) 

PAHs smaller than ~ 4.6 A will be photo destroyed in a time scale shorter than the age of the 
HD 141569A system (see §4.2). PAHs will also be efficiently removed from the disk by radiation 
pressure (see §4.3). Therefore, the PAH component in the HD 141569A disk must be continuously 
replenished, most likely through the evaporation of the icy mantles 9 of grains produced by the 
breakup of larger bodies. The evaporation of icy mantles takes place at r < 100 AU for micron- 
sized dust and at r < 30 AU for mm-sized dust (see Figure 6). Therefore, it is reasonable to assume 
that the free-flying PAH molecules, responsible for the observed PAH IR emission bands, are mainly 
distributed in the inner r < 100 AU region. If we assume that the spatial distribution of PAHs 
follows that of the porous dust at r < 100 AU, the mean starlight intensity to which the PAHs in 
the HD141569A disk are exposed is (U) » 1.24 x 10 6 (see Eq.[16]). 



3.1.4. PAH Ionization 

PAHs acquire charge through photoelectric emission and collisions with electrons and ions. 
Ideally, one can calculate the steady-state charge distribution for a PAH molecule of radius a from 
the balance between the electron capture rates and the photoelectron emission rates plus the ion 
capture rates (e.g. see Weingartner & Draine 2001). Since little is known regarding the electron 



(m* ~ 2.3m0; van den Ancker et al. 1998) and HD 141569B (m* ~ O.5m0; Weinberger et al. 2000) system with a mass 
ratio of q w 4.6, the Roche lobe radius is R CI « 0.49 A q 2/a / [o.6 q 2/3 + In (l + q 1/3 )] w 0.51 A, where A (« 990 AU 
[Weinberger et al. 2000]; « 1065 AU [Boccaletti et al. 2003]) is the physical separation between HD 141569A and 
HD 141569B. However, if HD 141569B is out of the HD 141569A disk plane, the physical separation A must be much 
larger than 990 AU. As will be discussed in §4.1, models with a larger r out value also provide good fits to the observed 
SED (see Figure 7). 

9 PAHs are abundant and widespread in the interstellar medium (see Leger & Puget 1984, Allamandola, Tielens, 
& Barker 1985, Li & Draine 2001a). It is expected that the icy mantles of the grains in dust disks, formed from 
condensation of volatile molecular species in the dense cloud phase, contain a substantial amount of PAHs since 
the freeze-out of PAHs onto the icy mantles are also expected to occur in the dense cloud phase and during the 
aggregation process. 
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density n e and its distribution in the HD 141569A disk, we will not carry out detailed calculations of 
the PAH charge distribution; instead, we will just calculate and compare the photoionization rates 
(ki on ; see §A) and the electron recombination rates (k rec ; see §A and §B) for PAHs as a function of 
size and as a function of distance from the central star. 

In Figure 2 we present the photoionization time scales T loIi (= l/k- loIi ) for PAHs of various sizes 
at 3 characteristic distances r = 20, 55, 90 AU (corresponding to the maximum and half-maximums 
of the PAH distribution; see Eq.[13]), calculated from the method described in §A. It is seen in 
Figure 2 that T\ on decreases as PAHs become larger and as PAHs are closer to the central star. The 
former is because when PAHs become larger, their ionization thresholds Ejp decrease (see Eq.[A6]) 
and their UV/visible absorption cross sections increase so that their photoionization rates increase 
(see Eq.[A7]). The latter is because, closer to the illuminating star, the UV radiation field is more 
intense (r; on oc r 2 ; see Eq.[A7]). 

Assuming that the source of the electrons in the HD 141569A disk is dominated by the cosmic 
ray ionization of H2, and adopting a total H2 mass of mn 2 w 115 m® inferred from submillimeter 
measurements of CO (Zuckerman, Forveille, & Kastner 1995), we calculate the electron recombina- 
tion time scales r rcc (= l/k ICC ; see §B for details). As shown in Figure 2, the electron recombination 
time scales are by ~ 3 orders of magnitude smaller than r; on . Therefore, PAHs in the HD 141569A 
disk will be negatively charged. 

3.1.5. Model Parameters 
We have two sets of parameters to be specified or constrained: 

1. The parameters jtipah, 4>ion, oo 5 and o determining the abundance, ionization fraction and 
size distribution of the PAH "log-normal" population (see Eqs. [2,3,10]). 

2. The parameters a p , P, and a determining the abundance, the fluffiness and the size distri- 
bution of the porous dust population (see Eqs. [6, 9, 10, 13, 14,15]). 

For a given set of ao and a, the total amount of PAH dust jtipah is determined by the absolute 
level of the observed PAH flux density (see Eqs. [2,3, 10]). Since little is known regarding the electron 
density n e and its distribution in the disk, it is impossible to perform accurate calculations on 
0ion- But we have seen in §3.1.4 that PAHs in the HD 141569A disk will be negatively charged 
by capturing electrons from the cosmic ray ionization of H2. Since we do not distinguish the IR 
properties of PAH - anions from PAH + cations (Li & Draine 2001a), 10 it is appropriate to take 
4>ion = 1- But we will also consider models with 4>- lon = 0.5,0 (see §4.1 and Figures 8b, c). 



10 It has been shown that the IR properties of PAH anions closely resemble those of PAH cations (e.g., see Szczepan- 
ski, Wehlburg, & Vala 1995; Langhoff 1996; Hudgins et al. 2000) except for the very strong 3.3 /tm C-H stretch 
enhancement in the anion (Szczepanski et al. 1995; Hudgins et al. 2000). However, Bauschlicher & Bakes (2000) 
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As shown in Li & Greenberg (1998) and Li & Limine (2003), cometary-type porous dust with 
P ~ 0.90 is able to reproduce the observed SEDs from the IR to millimeter for the disks around 
the MS star (3 Pictoris and the pre-MS star HR4796A; therefore, we will start with P = 0.90 in 
modeling the HD 141569A SED; but models with larger P or smaller P will also be discussed (see 
§4.1 and Figure 7). 

It is interesting to note that, as early as 1965, based on a Monte Carlo simulation, Cameron 
Sz Schneck (1965) showed that a very open structure ("fairy castle") with a porosity in the range 
of 0.83 < P < 0.89 (Cameron & Schneck [1965] used the term "underdense") is expected for 
dust aggregates assembled from randomly incident particles, depending on their shapes, sizes and 
incidence angles (also see Cameron 1995). A porosity of P « 0.90 for the "cold-coagulation" dust 
(which implies a mass density of ~ 1.7 x [1 — P'] gem" 3 ~ 0.46 gem -3 for this dust assuming ice 
mantles coated on the individual subgrains [P' « 0.73]; see Appendix B in Li & Lunine 2003) is 
consistent with the mean mass density « 0.5 gem -3 of cometary nuclei determined from studies of 
comet splitting and of non-gravitational forces (Rickman 2003; Whipple 1999). Comets are likely 
to have been formed initially through the cold accumulation of interstellar dust (e.g. see Greenberg 
& Li 1999). Hence, the ice-coated "cold-coagulation" dust aggregates we consider here are plausible 
building blocks of cometesimals and, hierarchically, comets. 

Similar to the PAH mass parameter mpAH, for a given set of P and a, the total amount of 
porous dust (characterized by cr p ; see Eqs. [11,13,14,15]) is determined by the absolute level of the 
observed flux density at wavelengths A > 10 jim. 

Therefore, we are left with N pax = 3 adjustable parameters: ao, cr, and a. 

For a model with N par adjustable parameters, the goodness of fit is measured by 

d.o.f . iV par + 105 

where [A-F\] mo d is the model spectrum (see Eq.[10]), and [AF\] bs is the (dereddened, stellar pho- 
tospheric contributions subtracted) observational spectrum which includes the 10 IRAS, OSCIR, 
MIRLIN, and SCUBA photometric data points and the 95 CGS3 spectroscopic data points (see 
Table 1). 



predict that PAH" anions have band strengths intermediate between those of neutrals (strong 3.3 /im C-H stretching 
and 11.3, 11.9, 12.7 fj,m out-of-plane C-H bending modes) and PAH + cations (strong 6.2, 7.7 /im C-C stretching and 
8.6 /im C-H in-plane bending mode). In any case, the detailed differences in the band strengths between PAH - anions 
and PAH + cations would not alter our conclusions since, as shown in §4.1, models containing a fraction (< 50%) of 
neutral PAHs are also capable of fitting the observed PAH mid-IR spectrum (see Figures 8b, c). 



(17) 
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3.2. Model Results 

We first consider the "cold-coagulation" dust model. As shown in Figure 3, our best-fit model 
(no. 1) provides an excellent fit to the entire spectrum (including both the broadband photometry 
and the PAH emission features) with emission from a mixture of PAH cations (0j on = 1) and 
porous dust with a porosity of P = 0.90, a power-law size distribution index of a « 3.3, and a 
total dust mass of ma ~ 3.56 m®. The PAH component is characterized by ao ~ 2.5 A, a ~ 0.3, 
and mpAH ~ 7.10 x 10~ 6 to©. See Table 2 for a full set of model parameters. For the sake of 
convenience, hereafter we call this model the "canonical cold-coagulation model". 

To illustrate the contributions to the observed IR emission from dust in different regions, we 
show in Figure 3a the IR emission from warm dust in the inner region of r < 105 AU, and from cool 
dust in the inner (105 < r < 235 AU) and outer rings (235 < r < 500 AU). While the cool dust in 
the two rings dominates the emission at A > 60 /Ltm and accounts for f» 80% of the 25 fim IRAS flux 
(it is the dust in the outer ring that is responsible for the 1350 fim SCUBA emission), the mid-IR 
10.8 and 18.2 ^m radiation detected by the OSCIR/Keckll instrument (Fisher et al. 2000; plotted 
as filled-triangles in Figure 3) are exclusively emitted by warm dust at r < 105 AU. This confirms 
earlier work of Fisher et al. (2000) and Augereau et al. (1999a) who suggested the existence of two 
separate populations of dust grains in the HD 141569A disk: inner, warmer grains that emit the 
mid-IR (10 and 18 /im) radiation, and more distant grains that are responsible for the scattered 
near-IR flux seen in the NICMOS image and the far-IR (60 and 100 /im) emission detected by 
IRAS. It is worth noting that in our model, the r < 105 AU warm dust fully accounts for both the 
10.8 [im and the 18.2 /jm emission, with very little contribution from the PAH component (also see 
Figure 3b). This is consistent with the discovery of Fisher et al. (2000) that the emitting regions 
at both wavelengths are of comparable size and therefore the same population of dust may emit 
both the 10.8 and 18.2 /im radiation. 

Our model yields a total dust IR flux of f°° ° F\d\ ps 5.89 x 10~ 10 erg s _1 cm" 2 , corresponding 

912 A 

to an IR luminosity of ~ 0.18 L©, ~ 8.1 x 10 of the total stellar luminosity (L* « 22.4 L©). The 
fractional contributions of the PAH component, the inner warm dust component, the inner ring, 
and the outer ring (see Figure 3a), are approximately 12%, 26%, 25%, and 37%, respectively. 

Our model predicts a vertical optical depth at visible wavelengths of t(V) ~ 0.024 at r p = 
203 AU (see Table 2). The vertical optical depths at other radial distances are smaller than this 
as can be seen from Eqs. (12-15). The in-plane optical depth is also much smaller than one. 11 
This justifies the optically-thin treatment employed in the entire paper. The PAH component 
causes little extinction: its contribution to the visual vertical optical depth at r = 55 AU where the 
distribution of PAHs in the HD 141569A disk peaks (see §3.1.3) is only r(V) « 1.54 x 10" 4 . 12 



n The in-plane optical depth can be calculated from t\ = / r r ° ut cr(r)/ [\/2ivH(r)] dr / a amax Cg^(a, A) (dnpo/da) da 
where H is the vertical scale height (see Eq.[B4]). Our best-fit model leads to r(V) w 0.097. 

12 At a first glance, one may argue that the PAH optical depth at visual wavelengths would also be about 12% of 
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Similar results are obtained for the "hot-nebula" dust model. We show in Figure 4 the best- 
fit P = 0.90 (P' ~ 0.80) "hot-nebula" dust model spectrum (model no. 2; see Table 2 for model 
parameters). In comparison with the observed SED, the overall fit is acceptable except that the 
model predicts a strong crystalline silicate emission band at 21 /im which is not seen in the CGS3 
spectrum (Sylvester et al. 1996), suggesting that only a small fraction of the hot silicate dust is in 
crystalline form. 

One should keep in mind that the "cold-coagulation" model and the "hot-nebula" model rep- 
resent two extremes; the actual dust in the HD 141569A disk is likely to be somewhere intermediate 
between these two types. Therefore, we consider a dust model consisting of a mixture of these two 
dust types with various mass mixing ratios. 13 Let /hn be the mass fraction of the "hot-nebula" dust 
in the mixture of highly-processed dust and "unaltered" dust. Let (XF\/m,i) cc and (XFx/rrid)^ re- 
spectively be the IR flux for per gram "cold-coagulation" dust (model no. 1) and "hot-nebula" dust 
(model no. 2). Let m^ c and m^ n respectively be the dust mass required for the "cold-coagulation" 
model (no. 1) and the "hot-nebula" model (no. 2) to reproduce the observed SED. The model spec- 
trum expected from a mixture of these two kinds of dust with a mixing mass fraction of / nn is 
simply 



(l-/hn)(^) +fJ K! ' 



(l-/h„)m^ c + / hn m^ . (18) 



m d/cc V m d J hnJ 

In Figure 5 we plot the model spectra with /hn = 10%, 20% and 30%. As can be seen in Figure 
5b, the non-detection of the 21 fim crystalline silicate feature places an upper limit of /h n ~ 10%. 

Finally, we present in Figure 6 the equilibrium temperatures for the best-fit "cold-coagulation" 
dust (P = 0.90, P' « 0.73) and "hot-nebula" dust (P = 0.90, P' « 0.80) as a function of grain size 
at a distance of r = 55, 203 AU from HD 141569A, and r = 70 AU from HR4796A. 



4. Discussion 

The narrow, ring- like structure was originally suggested by Sylvester & Skinner (1996) for 
the HD141569A disk in modeling the IRAS data (plus the JCMT upper limit at A = 1.1mm). 
However, their best-fitting models required a huge inner cavity with a radius of > 670 AU. On 



the total since the PAH component emits ~ 12% of the total flux so that it must have absorbed ~ 12% of the total 
stellar radiation. We remind the reader that in our model the PAH dust lies closer to the star (r < 105 AU) so that 
it is exposed to a radiation field w 1.24 x 10 6 /4.03 x 10 4 « 31 times stronger than that of the porous dust component 
(see Eq.[16] and §3.1.3). Another reason lies in the fact that PAHs absorb more efficiently in the UV compared to 
the visible than do micron-sized silicate or carbon grains (see Figure 2 in Li & Draine 2001a, Figure 1 in Li & Draine 
2002a, Figures 4a and 5a in Draine & Lee 1984). 

1,5 For simplicity, we assume that they are separate populations. In reality, it is more likely that the highly processed 
"hot-nebula" dust and the "unaltered" protostellar interstellar dust are mixed to some extent and hence form porous 
heterogeneous aggregates. 
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Table 2: Models for the HD 141569A dust disk IR emission* 



model 


dust 


P 


p ,b 




M d 


(«=»>«• 




m d 




mpAH* 


t p (V)S 


X 2 /Af 


note 


no. 


typo 








(pm) 


(Mm 2 ) 


(10~ 8 cm 3 ) 


(10 28 g) 


(10 4 cm" 2 ) 


(10 22 g) 


(io- 2 ) 






1 


cold-coag. 


0.90 


0.73 


3.3 


1.77 


7.18 


0.21 


2.13 


4.84 


4.24 


2.41 


1.89 


preferred 


2 


hot-ncbu. 


0.90 


0.80 


3.0 


2.00 


18.4 


2.00 


6.63 


1.52 


5.07 


1.85 


3.67 




3 


cold-coag. 


0.90 


0.73 


3.2 


0.18 


0.099 


0.0027 


2.52 


430.5 


4.29 


2.18 


2.27 


a-max —0.1 fJ.m 


4 


cold-coag. 


0.90 


0.73 


3.6 


16.2 


426.5 


9.65 


1.66 


0.081 


5.51 


2.21 


1.72 


Umax = 10 M m 


5 


cold-coag. 


0.90 


0.73 


3.3 


1.77 


7.18 


0.21 


3.01 


4.66 


4.32 


2.32 


1.70 


rout = 1000 AU 


6 


cold-coag. 


0.95 


0.87 


3.2 


1.83 


9.26 


0.44 


1.73 


3.89 


4.28 


2.38 


2.00 




7 


cold-coag. 


0.80 


0.46 


3.5 


1.67 


4.95 


0.049 


1.38 


6.58 


4.22 


2.21 


2.13 




8 


cold-coag. 


0.90 


0.73 


3.3 


1.77 


7.18 


0.21 


2.13 


4.84 


4.45 


2.41 


2.34 


a = 3.5 A, a = 0.4 


9 


cold-coag. 


0.90 


0.73 


3.3 


1.77 


7.18 


0.21 


2.13 


4.84 


6.19 


2.41 


4.97 


an — 3 A, a — 0.2, ion = 0.: 


10 


cold-coag. 


0.90 


0.73 


3.3 


1.77 


7.18 


0.21 


2.13 


4.84 


9.95 


2.41 


24.4 


an = 3 A, ct = 0.2, ion = 


11 


cold-coag. 


0.90 


0.73 


3.3 


1.77 


7.18 


0.21 


2.13 


4.84 


4.47 


2.41 


2.02 


«mt„ H =4-6 A 



"Unless otherwise stated, all models assume (1) a m in = 1 pm and a ma x = 1cm for the porous dust component; (2) 
a = 2.5 A, a = 0.3, a m f" = 3.5 A, and ion = 1 for the PAH component; (3) r G [0.2, 500 AU] for the dust spatial 
distribution given in Eqs. (13-15). 

6 The porosity of the ice-coated porous aggregate (see Appendix B in Li & Lunine [2003]). 

c The index of the power-law size distribution for the porous dust component: dn(a)/da ~ a~ a for a G [a m ; n , a mra ] 
(see Eq.[6]). 

d (a 7 ) = J^* max da a? dn(a) /da/ J^ max da dn(a)/da. The dust surface mass density can be written as a(r; m) = (m) a(r) 
with the mean dust mass (m) = (47r/3)(a 3 )(p) where (p) « 2.5(1 - P)gcm~ 3 , (p 1 ) « 1.7(1 - P')gcm~ 3 for the 
"cold-coagulation" model; and (p) « 3.5(1 - P)gcm~ 3 , (p 1 ) w 2.4(1 - P')gcm -3 for the "hot-nebula" model (see 
Appendix B in Li & Lunine [2003]). 

e The mid- plane surface density (z — 0) at r p = 203 AU for the porous dust component (see Eq.[14]). 

5 We assume that the spatial distribution of the PAH component follows that of the inner warm porous dust (which is 

responsible for the 10.8 and 18.2 pm mid-IR emission observed by Fisher et al. [2000]; see Figure 3a and Eq.[13]) at 

r < 105 AU (see §3.1.3). 

9 The vertical optical depth at A = 0.55 pm and r p = 203 AU. The contribution to t(V) from the PAH component is 
negligible (see §3.2). 

the other hand, it was shown by Malfait, Bogaret, & Waelkens (1998) that a flat, conventional 
disk extending from ~ 7AU to ~ 1400 AU with a power-law index of —0.6 for the dust spatial 
distribution was also able to fit the IRAS data. 

Malfait et al. (1998) modelled the HD 141569A IRAS data in terms of dust with a hypo- 
thetical opacity law of K a b s oc A -1 and a hypothetical radial-dependent temperature of Tk(r) = 
6180 (r/i?*) -2 / 5 K. Their inferred dust spatial distribution of dn(r)/dr oc r~ 0,6 with r G [7,1400AU] 
is inconsistent with the current view of a central void at r > 20 AU (e.g. Marsh et al. 2002; Mouillet 
et al. 2001). In the work of Sylvester & Skinner (1996), the dust in the HD 141569A disk was mod- 
elled either as a mixture of two separate components (amorphous silicate and amorphous carbon) 
or as a single pure amorphous carbon component. The dust was taken to be compact spheres with 
a power-law size distribution of dn/da ~ a -3 ' 5 in the size range of 50A < a < 100 t(m. The vari- 
ation of dust density within the disk with distance from the star was also treated as a power-law. 
Their model did not treat the PAH emission bands. It is therefore not surprising that their model 
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spectrum was deficient at A < 15 /mi. 

Using the OSCIR/Keckll 10.8 and 18.2 //m flux ratio, Fisher et al. (2000) placed an upper 
limit of a max ~ 1 /im on the radius of the dust grains responsible for the mid-IR emission at both 
wavelengths if they are made of compact silicate spheres. Since this conclusion was drawn from 
the assumption that the dust is located at a distance of 20 AU from the star, a wide range of grain 
sizes may be allowed if the mid-IR emitting dust is wide spread within r < 105 AU (see Eq.[13]). 

In contrast, a lower limit of a m j n ~ (0.6 ± 0.2) //m/ (1 — P) was derived by Boccaletti et al. 
(2003) for grains with a porosity of P, assuming that the slight surface brightness difference between 
their ground-based 2.2 /im imaging data and the HST NICMOS data at 1.1 fxm (Weinberger et al. 
1999) and 1.6 /zm (Augereau et al. 1999a) is caused by the color effect of the grains. A grain size 
distribution steeper than dn/da oc a~ 3 was also inferred from such an analysis. However, these 
conclusions are very sensitive to the scaling factor between the target star and the calibrator star 
which is uncertain. Once the scaling factor uncertainty is included in their analysis, a much wider 
range of grain sizes is allowed. 

In comparison with the narrow, sharply bounded HR4796A disk (Schneider et al. 1999), the 
dust annulus around HD 141569A is about 9-10 times wider and displays a more complex ring-gap- 
ring morphology (Augereau et al. 1999a; Weinberger et al. 1999; Clampin et al. 2003). In addition, 
for the HD 141569A disk, there is clear evidence for the existence of a warm dust population at 
r < 105 AU inward of the inner ring (Marsh et al. 2002). While a warm "zodiacal" dust component 
at a radial distance of a few AU from HR4796A was invoked by Koerner et al. (1998) and Augereau 
et al. (1999b) to account for the emission at A < 12 /mi, it was shown later by Li & Lunine (2003) 
that the porous dust model naturally explains the entire dust emission including that at A < 12 /im 
and therefore, there does not seem to be any strong evidences for a "zodiacal" dust component in 
the inner disk of HR4796A. Very recently, Marsh et al. (2002) found that the radial optical depth 
profile of the HD 141569A disk at A = 12.5 /im appears to increase inside of r = 20 AU at a 2a level. 
Does this indicate the existence of a "zodiacal" dust component in the HD 141569A disk? SIRTF 
imaging at the 3.6, 4.5, and 5.8 /jm IRAC bands (see Table 3 and §4.4) may be able to address this 
question since the "zodiacal" dust may also reveal its existence by emitting at the IRAC bands. It 
is interesting to note that the M band (4.77 //m) photometry does appear to show a small excess 
over the stellar photosphere (see Figure la). But the uncertainties both in the photometric flux 
determination and in the Kurucz model atmosphere representation prevent a definite conclusion. 

The time for complete sublimation of the ice mantles coated on the individual tenth-micron sub- 
grains of a porous aggregate of temperature T (see Figure 6) is r su bi ~ 1.5xl0 -12 (Aa//jm) (picc/g cm -3 ) 10^ 
(Backman & Paresce 1993) where Aa ~ 0.04 /im is the ice mantle thickness. By integrating 
r subi over the dust size distribution dn-pn/da, we estimate the mean sublimation time scales to be 
(^subi) ~ 7.6 x 10 5 , 3.3 x 10 32 , 2.8 x 10 10 yr respectively for the best-fit "cold-coagulation" dust at 
a distance of r = 55 AU, 203 AU from HD141569A, and 70 AU from HR4796A. If PAHs in dust 
disks indeed originate from the sublimation of icy mantles, we expect to see PAHs in the inner 
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warm region of the HD 141569A disk (r < 100 AU) where r su bi < r agc , but not in the ring regions 
of the HD 141569A (r > 100 AU) and HR4796A disks (r > 70 AU) where r subl > r agc . 

The 3.3, 6.2, 7.7, 8.6 and 11.3 /xm vibrational bands diagnostic of PAH molecules are not seen 
in the HR4796A disk. As mentioned above, this is mainly because in the r = 70 AU ring region, 
the sublimation of icy mantles which ejects PAHs is too slow (r su bi 3> r agc ~ 8 ± 3Myr). Even 
if PAHs are produced by mechanisms other than sublimation (e.g. shattering of carbonaceous 
grains), the recondensation of PAHs in the icy mantles of the porous dust component occurs 
on a time scale of a few thousand years, much shorter than the age of the HR4796A system: 
7-cond ~ (vpAH"- d vr(a 2 )) _1 ps (vr(mpAH)/8/cTg as ) 1/2 {^/^H/a^) {^{a 2 ))' 1 ps 1570yr where T gas (r) ps 
566 (r/AU) -1 / 2 K ps 67K is the gas temperature; (topah) ~ 7.99 x 10~ 22 g is the average mass of 
the PAH molecules with a size distribution like that of the HD 141569A disk (i.e., ao = 2.5 A and 
a = 0.3; see §3.2); H ps 0.5 AU is the vertical scale height (Kenyon et al. 1999); a p ps 4.90 x 10 4 cm" 2 
(Li & Lunine 2003) and na ~ <j p /V2itH are respectively the mid-plane surface density and volume 
density of the porous dust component at r p = 70 AU; vr(a 2 ) ps 9.02 x 10~ 7 cm 2 (Li h Lunine 
2003) is the total grain surface areas integrated over the dn/da oc a~ 2 - 9 size distribution with 
a £ [1 /xm, lcm] (see Table 1 [model no. 2] in Li & Lunine 2003). Similarly, the time scale for the 
accretion of PAHs on grains at r = 203 AU from the HD 141569A disk is about 9 x 10 5 yr, shorter 
than the age of the HD 141569A system. Therefore, free-flying PAH molecules are unlikely to exist 
in the ring regions of the HD 141569A and HR4796A disks. 

To summarize, the respective presence and absence of free-flying PAHs in the HD 141569A and 
HR4796A disks can be attributed to the structural differences between these two disks. It is worth 
noting that the (3 Pictoris disk, modelled best as a wedge-shaped disk with a constant opening 
angle (Artymowicz, Burrows, & Paresce 1989, Backman, Gillett, & Witteborn 1992), is more like 
the HD 141569A disk than the HR4796A disk does in the sense that the (3 Pictoris disk also extends 
to the inner warm region where the sublimation of icy mantles occurs. However, the existing mid-IR 
spectra of the j3 Pictoris disk do not appear to show any evidence of PAHs (Telesco & Knacke 1991, 
Knacke et al. 1993, Pantin, Waelkens, & Malfait 1999, Weinberger, Becklin, & Zuckerman 2003). 
It is not clear whether the absence of PAHs in the (3 Pictoris disk is due to its evolutionary status 
or just because the PAH emission features are swamped by the strong 9.7/xm and 11.3 /xm silicate 
bands. With an age of ps 12^ Myr (Barrado y Navascues et al. 1999, Zuckerman et al. 2001), the 
(3 Pictoris disk is more evolved than the HD 141569A disk. SIRTF spectroscopy of (3 Pictoris at 
A < 8/xm will allow us to draw a more definite conclusion regarding the presence or absence of 
PAHs in the inner disk of f3 Pictoris. 

The presence of transiently heated PAH molecules in the dusty environments of Herbig Ae /Be 
stars has been observationally well established (see e.g. Whittet et al. 1983; Brooke, Tokunaga, & 
Strom 1993; Sylvester et al. 1996; Siebenmorgen et al. 2000). Natta Sz Kriigel (1995) have calculated 
the mid-IR emission spectra expected for PAHs around Herbig Ae/Be stars, assuming a spherical 
shell geometry for the dust distribution. Sylvester, Skinner, & Barlow (1997) have modelled the 
SEDs of 8 Vega-like dust disks with temperature-fluctuating small silicate and amorphous carbon 
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grains of a = 3 A and 5 A included. So far, however, little has been done in modeling the PAH 
component in the HD 141569A disk. This work actually represents the first successful modeling of 
its entire SED, including the PAH emission features. 

4.1. Robustness 

The dust model presented here is fully described by 10 parameters: (1) aJ^ n H , ao, <J, and 0i on 
for the PAH component; (2) a m i n , a max , a, and P for the porous dust component; (3) r- m and r out 
for the dust spatial distribution dn{r)/dr. But, as discussed in §3.1.5, we actually have only 3 free 
parameters: ao, a, and a. 

Other parameters have been constrained by rather general considerations: (1) we take = 
3.5 A, the radius of the smallest survival interstellar PAH molecule (§3.1.1); (2) we take </>i on = 1 
since PAHs in the HD 141569A disk are expected to be negatively charged (§3.1.4); (3) we take 
a min = 1 A* m since this is the smallest size that makes sense in the framework of interstellar dust 
aggregation (see Footnote-5); (4) we take a max = 1cm; this size is not well constrained but neither 
is it crucial, since grains larger than ~ 100 /mi emit like blackbodies and their IR emission spectra 
are size-insensitive (see Figure 6); (5) we take P = 0.90 since dust of such a high porosity reproduces 
the SEDs of the [3 Pictoris disk and the HR4796A disk; a porosity of P « 0.90 is expected for 
fluffy aggregates formed by the accumulation of interstellar dust and is also consistent with the 
low density nature of cometary nuclei (§3.1.5); (6) we take n n to be the radial distance where the 
refractory dust starts to evaporate (§3.1.3); (7) we take r ou t = 500 AU which is expected from disk 
truncation caused by the tidal effects of HD 141569B (Footnote-8). The dust spatial distribution 
is taken to be that derived from the near-IR imaging of scattered starlight and mid-IR imaging of 
dust thermal emission (Eqs. [13-15] and §3.1.3). 

To be complete, we have also considered models with parameters differing from the nominal 
ones described above. In the following, except the specifically stated parameters, all other param- 
eters remain the same as those of the canonical "cold-coagulation" model derived in §3.2 (model 
no. 1; see Figure 3). See Table 2 for all the key parameters. 

• The a m i n = O.l^m model (no. 3; a = 3.2; Figure 7; red solid line) — is almost identical to 
the a m i n = 1 /im model (no. 1; Figure 3) except its slight deficiency at the IRAS 60 /mi band. 

• The a m i n = 10 /xm model (no. 4; a = 3.6; Figure 7; blue dotted line) — emits too much at the 
IRAS 60 and 100 /mi bands, and too little at the MIRLIN/Keckll 17.9, 20.8 /im bands and 
the IRAS 25 /mi band. 

• The r out = 1000 AU model (no. 5; a = 3.3; Figure 7; magenta dashed line) — is almost 
identical to the r out = 500 AU model (no. 1; Figure 3) since the cool dust at 500 < r < 1000 AU 
emits very little in the wavelength range of interest here. It emits only ~ 10% more at 
A > 200 /im than the r out = 500 AU model. 
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• The P = 0.95 model (no. 6; a = 3.2; Figure 7; green long dashed line) — is almost identical 
to the P = 0.90 model (no. 1; a = 3.3; Figure 3) except its slight deficiency at the SCUBA 
1350 /mi band. Its emission at 1350 /mi is ~ 50% lower than that of the P = 0.90 model. 

• The P = 0.80 model (no. 7; a = 3.5; Figure 7; cyan dot-dashed line) — is almost identical to 
the P = 0.90 model (no. 1; a = 3.3; Figure 3) except that it emits ~ 20% less at 100 < A < 
600 ^m. 

In summary, although the canonical "cold-coagulation" dust model presented in §3.2 is pre- 
ferred, some flexibilities are allowed in modeling the HD 141569A SED; i.e., no fine tuning is needed 
in order for our models to achieve an excellent fit to the observed entire SED. Therefore, the porous 
dust model made of interstellar materials, previously shown successful in reproducing the SEDs of 
the dust disks around f3 Pictoris and HR4796A, is also robust in modeling the IR emission from 
the HD 141569A disk. 

We have also tried to fit the observed SED using the interstellar PAH size distribution, char- 
acterized by ao = 3.5 A and a = 0.4 (Li & Draine 2001a). As shown in Figure 8a, the model 
spectrum predicted from the interstellar mixture of PAH cations (model no. 8) is also in an ex- 
cellent agreement with the observed mid-IR spectrum. In comparison with the canonical best-fit 
mixture (ao = 2.5 A and a = 0.3; model no. 1; see Figure 3), the interstellar mixture, rich in rela- 
tively large PAHs (see Figure 9), emits a little bit more (less) at wavelengths longward (shortward) 
of ~ 10 /im. This is because large PAHs do not reach temperatures as high as do small PAHs when 
heated by energetic photons. 

So far, all models assume (pi on = 1. We now consider models containing neutral PAHs. For 
illustration, two models are considered: a model with a mixture of both neutral and ionized PAHs 
(^ion = 0.5; model no. 9), and a model only consisting of pure neutral PAHs (4>i n = 0; model 
no. 10). In Figure 8b and 8c we respectively show the best-fit model spectra calculated from the 
0ion = 0.5 model and from the neutral PAH model. The PAH component of both dust models has 
a size distribution of ao ~ 3.0 A and a ~ 0.2 (see Figure 9). Again, other parameters (see Table 2) 
are the same as those of the canonical "cold-coagulation" model (no. 1). An inspection of Figure 
8c would immediately lead us to rule out the neutral PAH model since this model produces too 
strong an 11.3 /im feature (as well as a 12.7/tm feature) to reconcile with the CGS3 spectrum of 
Sylvester et al. (1996). The fit to the CGS3 spectrum by the <pi on = 0.5 model is acceptable except 
that it also appears to emit a bit too much in the 11.3 /tm band, although the quality of the CGS3 
spectrum is inadequate to disprove this model. Note that the model-predicted 11.3 /tin feature is 
much sharper than the observed one since we adopt the mean interstellar PAH bandwidths (see 
Table 7 of Li h Draine 2001a). As a matter of fact, the integrated flux in the calculated 11.3 /tm 
feature is comparable to the observed. The determination of the band strength of the 11.3 /tm 
feature is also complicated by its adjacent continuum. In our dust model, the warm porous dust in 
the inner region (r < 105 AU) dominates the continuum underneath the 11.3 /im feature. Future 
observational studies (e.g. SIRTF) will provide more information on 0j on . But at this moment, 
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it is safe to state that a large fraction of the PAH molecules responsible for the observed mid-IR 
spectrum in the HD 141569A disk must be ionized. 



4.2. PAH Destruction 

Upon absorption of an energetic photon, small PAHs, with an insufficient number of internal 
vibrational modes in which to distribute this photon energy, may be dissociated by ejection of a 
hydrogen atom, a hydrogen molecule, and/or an acetylene molecule (C2H2). See §A for a detailed 
discussion of this photoprocess. Following the method described in §A, we have calculated the 
photodestruction rates (fcdcs) for small PAHs exposed to stellar UV photons in the HD 141569A 
disk, represented by the photoejection rates of an acetylene molecule. In Figure 10 we show the 
PAH destruction time scales (rd es = 1/fcdes r~ 2 ) for PAHs at r = 20, 55, 90 AU as a function of 
PAH size. Apparently, PAHs smaller than a ~ 4.6 A at r < 105 AU are photolytically unstable: 
they are expected to be photodestroyed in a timescale shorter than the lifetime of the HD 141569A 
system. Larger PAHs are more stable since they have a larger number of vibrational modes so 
that they can easily accommodate the absorbed photon energy and it is therefore less likely for 
photodissociation to occur. 

In order to maintain a stable distribution of small PAHs in the inner region (r < 105 AU) of 
the HD 141569A disk, there must exist a source continuously replenishing the small PAHs at a rate 
of 

/" 105AU r y> a I" A ^pah (47r/3)a 3 ppAH nq , 

"iPAH = / apAH(r)2irrdr / da— (19) 

Jr in A PAH d a r dcs (a,r) 

in mm 

where the PAH surface density distribution, ctpah(^), is taken to be the same as that of the inner 
warm porous dust (r < 105 AU; see §3.1.3 and Eq.[13]) 

ctpah (r) = °"pah exp { -4 In 2 [(r/AU - 55) /70] 2 } (20) 

where (Jp AH is the PAH mid-plane surface density at r = 55 AU. For the canonical "cold-coagulation" 
dust model (no. 1; §3.2; Figure 3) with a = 2.5 A, a = 0.3, a^ 1 = 3.5 A, and cip AH pa 1.17 x 
10 13 cm -2 , we arrive at mpAH ~ 1-11 x 10~ 3 m^/yr! However, the PAH photodestruction time 
scales Tdes derived above may have been substantially underestimated for two reasons: (1) the pho- 
toejection of an acetylene molecule is not equivalent to the complete destruction of a PAH molecule; 
(2) PAHs can be repaired by accreting carbon atoms and/or ions from the gas (e.g. see Allain et 
al. 1996); the presence of carbon ions in the HD 141569A disk was revealed by the detection of the 
[CII] 158 fim emission (see Lorenzetti et al. 2002). Therefore, it is very likely that the required 
PAH replenishment rate estimated above has been considerably overestimated. 

But even if we assume that the above derived PAH mass loss rate due to photodestruction is 
valid and there is insufficient replenishment implying that there will be no PAHs smaller than ~ 
4.6A , our model with aj^ 1 = 4.6 A is still able to reproduce the observed SED. This is demonstrated 
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in Figure 11 for the "cold-coagulation" dust model with aj^ 1 = 4.6 A (model no. 11). Therefore, 



the exact knowledge of the PAH lower-cutoff size is not critical in modeling the HD 141569A 
SED. 



4.3. Radiation Pressure and Poynting-Robertson Drag 

In addition to the gravitational attraction from the central star, grains in the HD 141569A disk 
are subject to (1) radiative outward propulsion as a consequence of momentum transfer from stellar 
photons, and (2) Poynting-Robertson drag which causes them to spiral toward the gravitational 
force center as a consequence of angular momentum loss through radiation (Burns, Lamy, & Soter 
1979; Backman & Paresce 1993; Kriigel 2003). 

The ratio of the radiative force to the gravitational force is 

a r s R * F a t C abs(a, A) + (1 - (g)) C sca (a, A)] dX 

Ma) = 4cGM,(4vr/3)a3p PAH < 21 > 

where G is the gravitational constant, M* (« 2.3m Q ) is the stellar mass, C sca (a, A) is the scattering 
cross section for a grain of size a at wavelength A, (g) ("asymmetry factor") is the average value of 
the cosine of the scattering angle. 

We have calculated /3rp as a function of dust size for the best-fit canonical "cold-coagulation" 
dust model (no. 1; see Figure 3 and §3.2). As shown in Figure 12a, for grains smaller than ~ 
10 — 20 ^m, the radiation pressure (RP) overcomes the gravitational attraction (i.e. /3rp > 1) 
and, therefore, these grains will be blown out from the HD 141569A disk. But we note that a 
stable distribution of small grains (< 10 — 20 ^m) in the disk is possible provided that they are 
continuously replenished by collisions of larger bodies and/or larger grains so that at any time the 
disk contains a substantial population of small particles (see Krivov, Mann, & Krivova 2000). 



r out , , Z" 10 ^ , rfnpp (47T/3) qV d 

w-RP = / a{r)2irrar / da — ; — — 

Jn n A min da r R p(r) 



The dust removal rate due to the radiation pressure (RP) expulsion can be estimated from 

da dn PD (4vr/3)ay D 
da 

™ ir) a 1900 (sJau) 3 ' 2 yr < 23 > 

where trp, the RP timescale, is assumed to be comparable to the local dynamical timescale Td yn = 
2ir/Q(r) where f2(r) = (GM^/r 3 ) 1 ^ 2 is the Keplerian frequency. For the best-fit canonical "cold- 
coagulation" model (no. 1; see Figure 3 and §3.2), the RP dust mass loss rate is 7.87x 10~ 6 m^/yr. 

We have also calculated the Poynting-Robertson (PR) drag timescale tpr(<2, r) (the time it 
takes for a grain of size a at a distance of r to fall into the central star) from 

/ x ° 2r2 (4tt/3) a 3 /9 PD f . 
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In Figure 12b we show that, for grains smaller than ~ 20 /mi at a radial distance of r p = 203 AU, the 
lifetime associated with the PR drag is shorter than the HD 141569A age and hence, these grains 
will fall into the central star even if they are not ejected by radiation pressure. For the inner warm 
dust component at r < 105 AU, grains in the range of 10 — 20 /im < a < 100 — 200 /im, although 
stable against radiation pressure ejection, will be removed from the disk by the PR drag. Since the 
best-fit "cold-coagulation" model (no. 1) has ~ 80% of the total surface areas in grains smaller than 
100 /im, these grains must be efficiently replenished by cascade collisions of planetesimals and larger 
grains. By integrating the PR dust removal rate (47r/3) a? ppr> / tpr(ci, r) over the whole size range 
and over the entire disk, we estimate the PR dust mass loss rate to be ~ 1.38 x 10 -8 m^/yr for the 
best-fit canonical "cold-coagulation" model (no. 1). Therefore, over the life span of HD141569A, 
roughly 39 m® of dust is lost by radiation pressure and Poynting-Robertson drag. 

For the PAH component, the ratio of radiation pressure to gravity /3rp 87.4 is almost 
independent of size. This is because PAHs are in the Rayleigh limit at UV/visible wavelengths 
where the absorption of stellar photons is such that C a b s (o, A) oc a 3 /A and (g) ~ (Bohren & 
Huffman 1983; Kriigel 2003). The PR time scale for PAHs at r = 55 AU is about 6000 yr, also 
independent of size for the same reason. The PAH mass loss rates due to radiation pressure and 
Poynting-Robertson drag are ~ 3.09 x 10" 8 m®/yr and ~ 1.95 x lO~ 9 m0/yr, respectively. This 
implies that a total amount of ~ 0.16 of PAH molecules needs to be replenished over the life 
span of HD 141569A. 

4.4. Predictions for SIRTF 

The Space Infrared Telescope Facility (SIRTF) will be capable of sensitive imaging using the 
Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 /im, and using the Multiband Imaging 
Photometer (MIPS) at 24, 70, and 160 jum. In Table 3 we show the band-averaged intensities for 
our preferred dust models. 

SIRTF will also be able to perform low-resolution 5-40 /im and high-resolution 10-37 /im spec- 
troscopic observations using the Infrared Spectrograph (IRS) instrument. IR spectroscopy and imag- 
ing will provide powerful constraints on the HD 141569A dust spatial distribution and its chemical 
composition. 

Our best-fit canonical "cold-coagulation" dust model (no. 1; Figure 3) predicts strong PAH 
C-C stretching features at 6.6 /mi and 7.7 /im, a weak in-plane C-H bending feature at 8.6 /im, a 
weak out-of-plane C-H bending feature at 11.3 /im, a very weak C-H stretching feature at 3.3 /im, 
very weak out-of-plane C-H bending features at 11.9 /im and 12.7/mi, and very weak C-C bending 
features at 16.4, 18.3, 21.2 and 23.1 /im (see Li & Draine 2001a for details). With neutral PAHs 
included (e.g., see Figure 8b for the <f>\ on = 0.5 model [no. 9] which is not ruled out in the present 
study), stronger 3.3 /im and 11.3 /im features are expected. SIRTF spectroscopy at A < 8 /im will 
provide an important test of the PAH model with respect to the PAH size distribution, charging, 
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and survival in the HD 141569A disk. 

The "cold-coagulation" dust model (no. 1) produces a broad, smooth amorphous silicate feature 
at 9.7 fim which forms a plateau underneath the 11.3 /^m PAH feature (see Figure 3 and §3.2). In 
contrast, the "hot-nebula" dust model (no. 2; see Figure 4 and §3.2) predicts two crystalline silicate 
features at 11.3 and 23 ^m. While the 11.3 /im feature is swamped by the PAH C-H out-of-plane 
bending mode which is also at 11.3 iim, the 21 iim feature is prominent and the CGS3 mid-IR 
spectrum readily places an upper limit of about 10% on crystalline silicate dust (see Figure 5 and 
§3.2). With high resolution spectroscopy, SIRTF will allow us to perform more detailed studies of 
the silicate emission features so as to infer the degree of processing which the HD 141569A dust has 
experienced. 

Table 3: Dust IR emission (Jy) integrated over SIRTF bands predicted for our preferred models. 



model 


IRAC 


IRAC 


IRAC 


IRAC 


MIPS 


MIPS 


MIPS 


no. 


3.6 fj,m 


4.5 iim 


5.8 /im 


8.0 /xm 


24 /xm 


70 /xm 


160 /Ltm 


1 


0.011 


0.018 


0.12 


0.41 


1.86 


4.45 


2.27 


8 


0.0098 


0.015 


0.090 


0.38 


1.94 


4.57 


2.37 


9 


0.034 


0.023 


0.14 


0.43 


1.91 


4.52 


2.32 


11 


0.0098 


0.016 


0.10 


0.40 


1.91 


4.55 


2.35 



5. Conclusion 

We have modelled the mid-infrared to submillimeter spectral energy distribution of the double- 
ring disk around HD 141569A in terms of the porous dust model previously shown successful in 
reproducing the SEDs observed for the disks around (3 Pictoris (Li & Greenberg 1998) and HR4796A 
(Li & Lunine 2003). The dust is modelled as porous aggregates of either unaltered or heavily 
processed interstellar materials. While the former ( "cold-coagulation" dust) is made of amorphous 
silicate dust, carbonaceous dust, H20-dominated ices, and vacuum, the latter ("hot-nebula" dust) 
consists of a mixture of crystalline silicates, ices, and vacuum (§3.1). For both dust types, a 
vacuum volume fraction of 90% is assumed (§3.1.5). The dust spatial distribution is taken to be 
that derived from imaging observations of scattered light and dust thermal emission (Eqs. [13-15]; 
§3.1.3). We take a simple power-law for the dust size distribution (dn/da oc a~ a ) in the size range 
of 1 itm < a < 1 cm. In addition, a population of PAH molecules with a log-normal size distribution 
characterized by a peak-parameter ao and a width-parameter of a is invoked to account for the 
observed 7.7 i/m and 11.3 /an "unidentified infrared" emission features that are generally attributed 
to PAHs. Our principal results are: 

1. The "cold-coagulation" dust model with a porosity of P = 0.90, a size-distribution power 
index of a « 3.3, and a total dust mass of f» 3.56 me, together with a population of 
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charged PAHs with ao ~ 2.5 A, a « 0.3 and a total mass of topah ~ 7.10 x 10 6 m e provides 
an excellent fit to the entire SED, including the PAH mid-IR emission features (§3.2; Figure 
3). 

2. Although the P = 0.90 "hot-nebula" dust model provides an overall good fit to the observed 
SED, it predicts a strong crystalline silicate emission feature at 21 fim which is not seen in 
the currently available observational spectrum (§3.2; Figure 4). This places an upper limit of 
~ 10% on the fraction of "hot-nebula" dust (i.e. crystalline silicates; §3.2; Figure 5). 

3. The PAH molecules in the HD 141569A disk are expected to be negatively charged (Figure 
2; §3.1.4; §A), assuming that the cosmic ray ionization of H2 is the dominant source of 
electrons (§B). But it has also been shown that models consisting of a mixture of both 
neutral and charged PAHs are also capable of reproducing the observed PAH IR emission 
features, provided that the fraction of neutral PAHs is smaller than 50% (Figure 8; §4.1). 

4. Approximating the photodestruction rate of a PAH molecule by the photoejection rate of 
an acetylene molecule, PAHs smaller than ~ 4.6 A will be photodestroyed in a time scale 
shorter than the age of the HD 141569A system (Figure 10; §4.2; §A). In order to have a 
stable distribution of small PAHs (3.5 < a < 4.6 A) in the HD 141569A disk, we require a 
continuous replenishment of this component at a rate of ~ 1.11 x 10 _3 m^/yr. Such a high 
replenishment rate may not be essential, since (1) the photoejection of a C2H2 molecule does 
not necessarily lead to the complete destruction of a PAH molecule, and (2) the repair of a 
PAH molecule with a broken C-C bond through reacting with gas-phase C atoms and/ions 
is ignored in deriving the PAH destruction time scale. Furthermore, we have shown in §4.2 
(Figure 11) that a model exclusively consisting of PAHs larger than ~ 4.6 A also achieves a 
good fit to the observed SED. 

5. Grains smaller than ~ 10 — 20 fim will be radiatively expelled from the disk; grains at r = 
55 AU in the size range of ~ 10 — 20 /im < a < 100 — 200 fim will also be removed from the disk 
due to the Poynting-Robertson inward spiralling drag (at a closer [further] distance from the 
star, larger [smaller] grains will be removed). Collisions of planetesimals must continuously 
replenish the dust in the disk at a rate of « 7.88 x 10~ 6 m^/yi (§4.3). PAHs will also be rapidly 
removed from the disk by radiation pressure and, to a less degree, by Poynting-Robertson 
drag. A replenishment rate of « 3.28 x 10~ 8 m ffi /yr is needed to maintain the PAH disk 
(§4.3). 

6. Spectroscopic and broadband photometric predictions are made for SIRTF observations (§4.4). 

We thank M.J. Barlow for providing us with the HD 141569A CGS3 mid-IR spectrum. We 
thank D.E. Backman, A.G.W. Cameron, C. Kulesa, V. La Page, S. Leach, T.W. Rettig, and M.D. 
Silverstone for helpful discussions, suggestions, and/or clarifications. A. Li thanks the University of 
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A. Photophysics of Polycyclic Aromatic Hydrocarbon Molecules in the 

HD 141569A Disk 

Following the absorption of an energetic photon, a PAH molecule has 3 major competing decay 
channels to relax its energy: emission; 14 ionization; and photodissociation. The photodissociation 
process, critical for small PAHs, has 3 major loss channels: the loss of a hydrogen atom, a hydrogen 
molecule, and an acetylene molecule (C2H2). 

We use the Rice, Ramsperger, and Kassel (RKK) theory (Forst 1973) for the photodissociation 
rate for a PAH molecule of internal energy E after excitation 

/ px\3Af-7 

k x (E) = , E>El; 

= 0, E<E%; (Al) 

where ./V is the number of atoms in the molecule; E£ is the activation energy for the loss of the 
species-x (= H, H2,C2H2) related to bond dissociation (i.e., the minimum energy required for the 
dissociation to occur); is the frequency factor for the x-loss channel which reflects the efficiency 
of intramolecular redistribution of the energy over the vibrational levels of the molecule's ground 
electronic state after internal conversion. In Table 4 we tabulate the E^ and v& parameters for the 
loss of H, H2 and C2H2 experimentally obtained by Jochims et al. (1994) for small PAHs. 

Table 4: Activation energy E& and frequency factor v& (Jochims et al. 1994) . 



loss channel 


El 




X 


(eV) 




H 


2.8 


1 x 10 16 


H 2 


2.9 


1 x 10 16 


C2H2 


2.9 


1 x 10 15 



The IR photon emission rate for a PAH molecule of internal energy E is 

where the vibrational temperature T is obtained from the vibrational energy E using the Debye- 
model approximation for the PAH enthalpy (see Eq.[15] in Li & Draine 2001). 



14 The omission process is dominated by fluorescence (transitions between vibrational states of same multiplicity) 
in the IR and part in the visible. Phosphorescence (transitions between vibrational states of different multiplicity) is 
less important. In this work we will thus only consider IR emission via fluorescence. 
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Assuming that IR emission, ionization, and dissociation are the only relaxation processes for 
a PAH molecule excited by an energetic photon, the photodestruction rate (approximated by the 
ejection of an acetylene molecule) kdes, and the dehydrogenation rate kdhy for PAHs at a distance 
of r from the central star can be obtained from following equations 15 

kdcs[r) -\2r) J hc/E c^ k n + k R2 + k C2U2 + k m hc/X ' (A ' j 



2r J Jhc/Ef k R + k U2 + k C2 u 2 + ^IR h c/X 

Following Allain, Leach, & Sedlmayr (1996), we adopt an analytical formula for the ionization 
yield Y lon which was based on the experimental data of pyrene and coronene (Verstraete et al. 1990) 

li on = , E < E IP ; 

0.1 (E - Etp) 

= { -t- t ^ , Ejp < E < 8.9eV ; 

8.90 eV -Eip 

= 0.25 (E/eV) - 2.125 , 8.9<£<lleV; 

= 0.0596 (E/eV) - 0.03 , 11 < E < 13.6 eV ; (A5) 

where E and Ejp (ionization threshold) are in units of eV. We adopt a simple expression for the 
ionization threshold (see Eq.[2] in Weingartner Sz Draine 2001 for Z = 0), 

where the PAH radius a is in unit of angstrom. The photoionization rate for a PAH molecule 
located at a distance of r from the central star is 

^"(t) t*^^- (A7 » 

The electronic recombination rate for a PAH cation is 

(A8) 

1 /2 

where n e is the electron density; v e = (8kT gas /Trm e ) 1 is the electron thermal velocity; m e is 
the electron mass; e is the electron charge; ip = e/a is the PAH electrostatic potential; k is 
the Boltzmann constant; T gas is the gas temperature. We will take T gas (r) = (R k /2r) 1 ^ 2 T c g m 
605 (r/AU) -1 ^ 2 K for the HD 141569A disk. In Eq.(A8) we have assumed that the sticking proba- 
bility for recombination of an electron on the PAH cation equals one. The term 2etp/mVg accounts 
for the "Coulomb focussing" effect (Spitzer 1941). 



akU 



15 For simplicity, we assume that the absorbed photon will either ionize the neutral PAH molecule or participate in 
the photodissociation and vibrational relaxation. 
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B. Electron Density in the HD 141569A Disk 

Molecular CO gas in the HD 141569A disk has been detected in emission through its rotational 
transitions at submillimeter wavelengths (Zuckerman, Forveille, & Kastner 1995) and through its 
ro-vibrational transitions at near-IR (Brittain & Rettig 2002; Brittain et al. 2003). Brittain & 
Rettig (2002) also reported the detection of emission in the HD 141569A disk. In principle, 
one can estimate the electron density from the H^j~ and CO measurements based on the following 
assumptions - (1) H^~ is produced mainly through the cosmic ray (CR) ionization of H2 to H^: 
H2 + CR — > H^ +e~ and the ion-neutral reaction: H^ + H2 — > H^ + H; (2) H^ is destroyed primarily 
by its recombination with an electron: Hj~ + e~ — > H2 + H and by the ion-neutral reaction with CO: 
H^ + CO — > H2CO + + H (McCall et al. 1998). Let n e , nu 2 , nco, and n R + be the number densities 
of electrons, H2, CO, and H3", respectively; ?cr be the cosmic ionization rate; k c and kco be the 
rate constants for the H^-electron recombination and the Kj~-CO reaction, respectively. One can 
derive the electron density from 

?CR^H 2 - k C oncon R + 

He = 1 

n 3 

Brittain & Rettig (2002) estimated the column densities of CO and N H + respectively to be 
Nco ~ 10 11 cm" 2 and -/V H + ~ 7 x 10 10 cm -2 from the CO and H3" emission spectra obtained 
for the HD 141569A disk with the NASA Infrared Telescope (IRTF) using the Cryogenic Echelle 
Spectrograph (CSHELL). 16 With ?CR = 3 x 10- 17 s _1 , k co = 2.0 x lO^cn^s -1 , k e = 1.8 x 
10~ 7 cm 3 s _1 (see McCall et al. 1998 and references therein), and the nominal Nn 2 /Nco ~ 10 4 (Thi 
et al. 2001), one finds n e < from Eq.(Bl)! This suggests that CO and H^ may be distributed 
in different regions. Brittain & Rettig (2002) suggest that Hjj~ may originate from the extended 
envelope of a gas giant protoplanet at ~ 7AU. 

Therefore, it is not appropriate to derive n e from the H^!~ data; instead, we will place an upper 
limit on n e from nn 2 assuming that the cosmic ray ionization of H2 is the major source of electrons: 
n c « T^cR^Ha where 5Myr is the stellar age. 17 Below we will discuss how the number density 
of molecular hydrogen «h 2 is obtained. 



16 This assumes that both CO and are contained in the 17 < r < 50 AU region around HD 141569A. 

17 Atomic gas has also been seen in the HD 141569A disk. A double-peaked Ha emission was detected by Andrillat, 
Jaschek, & Jaschek (1990) and Dunkin, Barlow, & Ryan (1997), indicating a rotating gas disk close to the star. Dunkin 
ct al. (1997) also detected circumstellar Nal and Call absorption lines. Recently, Lorenzetti et al. (2002) reported 
the detection of C + gas. However, little is known regarding the atomic gas column density and its contribution to 
the electrons in the disk. Since it is shown in §3.1.4 that the amount of electrons provided by H2 alone is sufficient to 
(negatively) charge the PAH molecules in the HD 141569A disk, inclusion of electrons from the cosmic ray ionization 
of atomic gas would further support our conclusion (i.e., a large fraction of the PAH molecules are [see §3.1.4 and 
Figure 2] and must be [see §4.1 and Figure 8] charged in the HD 141569A disk). Furthermore, an additional amount 
of electrons would be provided by the X-ray ionization of H2 . 



(Bl) 
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Assuming that the disk is in vertical hydrostatic equilibrium, the spatial distribution of molec- 
ular hydrogen is 

nw (r,z) = , exp 



2 \HJ 



(B2) 



+oo 



^(r) = / n U2 (r,z)dz (B3) 

J — oo 

where £Jh 2 is the surface density of H2 ; the vertical scale height H is 

g(r) = f Jjjgg V /2 r 3/2 (B4) 

where M*(ps 2.3m Q ) is the stellar mass, ran is the mass of a hydrogen atom, [i is the mean molecular 
weight (fi ps 2.34 for a mixture of H2 and He). The total molecular hydrogen mass in the r < 130 AU 
region around HD 141569A was estimated to be m^' ps 5 — 115 based on the IRAM (Institut de 
Radio-Astronomie Millimetrique) observations of the J = 2 — 1 and J = 1 — rotational transitions 
of CO (Zuckerman et al. 1995). 18 More recently, Brittain & Rettig (2002) found that the CO gas, 
pumped by strong stellar UV radiation, is confined to the 17 < r < 50 AU region with a total mass 
of « 10 19 g. Assuming the standard ratio of H2 to CO (about 10 4 ), this leads to a total H2 mass of 
ps 2 x 10~ 5 m®. 

The H2 mass inferred by Brittain h Rettig (2002) differs by a factor of ps 7 x 10 6 from that of 
Zuckerman, Forveille, & Kastner (1995). This is because the UV-pumped CO detected by Brittain 
& Rettig (2002) only represents material near the star/disk interface and is not indicative of the 
total gas mass. The UV radiation is quickly scattered by dust and may be influential for only a 
small annulus of the disk. Submillimeter measurements of CO can measure the cold regions inward 
of ~ 130 AU (Zuckerman et al. 1995), where most of the cold CO molecules are in the ground 
vibrational state to which the near-IR observations are not sensitive (Brittain et al. 2003; T.W. 
Rettig, private communication). Therefore, we will adopt the Zuckerman et al. (1995) results for 
the H2 mass. 

We assume that in the r < 130 AU region, the spatial distribution of H2 follows that of the 
dust (see Eq.[13]), i.e., 

cr H2 (r) = cr£ 2 exp j-41n2[(r/AU- 55)/70] 2 } , r in < r < 130 AU ; (B5) 

where cr^ 2 is the mid-plane surface density of H2 at r = 55 AU which can be derived from 
(ps 5 — 115 m®; Zuckerman et al. 1995) 

/•130AU 

= / 2m H o- H2 (r)27rrdr . (B6) 



18 Zuckerman et al. (1995) originally gave to-h' ~ 20 — 460 m®. But this was derived from the assumption that the 
distance from HD 141569A to the Earth is d ^ 200 pc. Adopting d ^ 99 pc (van den Ancker et al. 1998), mH 2 should 
be reduced by a factor of ^ 4. 
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Fig. 1. — (a): Spectral energy distribution of the HD 141569A system (see §2 and Table 1). Red solid line: 
the T c ff = 10000 K and \gg = 4.0 Kurucz model atmospheric spectrum; green open triangles: optical and near-IR 
photometry (Sylvester et al. 1996); cyan filled triangles: the 10.8 and 18.2 /im OSCIR/Keckll photometry (Fisher et al. 
2000); black filled squares: the 12, 25, 60 and 100 fim IRAS photometry (Sylvester et al. 1996); magenta open squares: 
the 12.5, 17.9, and 20.8 ^m MIRLIN/Keckll photometry (Marsh et al. 2002). (b): The CGS3 mid-IR spectrum (blue 
filled circles; Sylvester et al. 1996) as well as the OSCIR, MIRLIN, and IRAS photometry. All observational data in 
(a) and (b) have been dereddened, but the stellar photospheric contributions to the IR photometric and spectroscopic 
bands have not been subtracted, (c): The HD 141569A radiation field at a distance of r = 4.75 x 10 4 AU from the 
star (red solid line) which equals the 912 A - 1 fim general interstellar radiation field (blue dashed line; Mathis et al. 
1983; see §3.1.1). 
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Fig. 2. — Photoionization (ri on ; upper part of this figure; see §A) and electron recombination (r rcc ; lower part; see 
§A and §B) time scales for PAHs in the size ranges of 3 A < a < 20 A at r — 20 AU (dotted lines), 55 AU (solid lines) 
and 90 AU (dashed lines) . In calculating the electron recombination time scales, the electron density n c is estimated 
from the molecular hydrogen mass determined by Zuckerman et al. (1995), assuming that the cosmic ray ionization 
of H2 is the dominant contributor to the electrons in the HD 141569A disk. The upper axis indicates the number of 
C atoms iVc in a PAH molecule. 
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Fig. 3. — Comparison of the "cold-coagulation" dust model spectrum (red solid line; model no. 1) with the observed 
spectral energy distribution. The model consists of two dust components: (1) the PAH cation component (cj>i on = 1; 
blue dashed line) with a log-normal size distribution (and a lower-cutoff of a^J 1 = 3.5 A) characterized by a peak 
parameter of ao = 2.5 A, a width parameter of a — 0.3, and a total mass of ttipah ~ 7.10 x 10 -6 m®, and (2) the 
porous dust component (green dotted line) characterized by a porosity of P = 0.90 [P' ~ 0.73), a power- law 
index of a ~ 3.3 for the dust size distribution (a € [1/im, lcm]), a dust surface density of a p ~ 4.84 x 10 4 cm~ 2 , 
and a total mass of ma ~ 3.56 m©. We adopt the dust spatial distribution (0.2 < r < 500 AU) derived from imaging 
observations of scattered light and dust thermal emission (see Eqs. [13-15] and §3.1.3). Upper panel (a) also illustrates 
the regional contributions of the porous dust in the inner warm region (r < 105 AU; cyan long dashed line), the inner 
ring (105 < r < 235 AU; magenta dot-dashed line), and the outer ring (235 < r < 500 AU; black dot-long dashed 
line). The observational SED includes the 10.8 and 18.2 /^m OSCIR/Keckll photometry (cyan filled triangles; Fisher 
et al. 2000), the 12, 25, 60 and 100 /im IRAS photometry (black filled squares; Sylvester et al. 1996), the 12.5, 
17.9, and 20.8 /^m MIRLIN/Keckll photometry (magenta open squares; Marsh et al. 2002), and the CGS3 mid-IR 
spectrum (black filled circles in lower panel [b] which illustrates the PAH emission bands; Sylvester et al. 1996). All 
observational data in (a) and (b) have been dereddened, with the stellar photospheric contribution subtracted. 
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Fig. 4. — Same as Figure 3 but for the P = 0.90 (P' w 0.80) "hot-nebula" dust model (no. 2). This model predicts 
a broad prominent feature at 21 /im due to crystalline silicates which is not seen in the CGS3 spectrum (Sylvester et 
al. 1996). 
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Fig. 5. — Model spectra expected from a mixture of the "cold-coagulation" dust (model no. 1; Figure 3) with a 
mass fraction of (1 — /h n ) and the "hot-nebula" dust (model no. 2; Figure 4) with a mass fraction of fhn (red solid 
line - /h n = 10%; blue dotted line - f^ n = 20%; green dashed line - /h n — 30%). The non-detection of the 21 fim 
crystalline silicate feature limits the fraction of "hot-nebula" dust (i.e. crystalline silicates) to /h n w 10%. 
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Fig. 6. — Equilibrium temperatures for "cold-coagulation" dust of P — 0.90 (P w 0.73) at a distance of r = 55 AU 
(solid line), and r = 203 AU (thin solid line) from HD 141569A (upper panel [a]), r = 70 AU (solid line) from HR4796A 
(lower panel [b]), and for "hot-nebula" dust of P = 0.90 (P' « 0.80) at r = 55 AU (dashed line) and r = 203 AU 
(thin dashed line) from HD 141569A (upper panel [a]), and r = 70 AU (dashed line) from HR4796A (lower panel 
[b]). The horizontal dotted lines plot the black-body temperatures of Tbb ~ 81.6, 42.5, 67. 6K respectively at r = 55, 
203 AU from HD 141569 A (upper panel [a]), and r = 70 AU from HR4796A (lower panel [b]). The discontinuity 
at T ~ 120 K in the temperature profiles (as a function of grain size) results from ice sublimation: grains hotter 
than ~ 120 K do not maintain icy mantles and therefore attain higher temperatures than their icy counterpart since 
silicate and carbon materials have a much higher UV/visual absorptivity than H2 0-dominated ice. 
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Fig. 7. — Comparison of the observed SED with the theoretical spectra calculated from the "cold-coagulation" 
dust models with a m i n = 0.1 ^im (model no. 3; a = 3.2; red solid line), a m ln = 10 )im (model no. 4; a = 3.6; blue 
dotted line), r out — 1000 AU (model no. 5; a = 3.3; magenta dashed line), P — 0.95 (model no. 6; a = 3.2; green long 
dashed line), and P = 0.80 (model no. 7; a — 3.5; cyan dot-dashed line). Unless otherwise specified, the remaining 
parameters of these models (see Table 2) are the same as those of the canonical "cold-coagulation" model (no. 1) 
described in §3.2 (see Figure 3). 
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Fig. 8. — Best-fits to the CGS3 mid-IR spectrum provided by (a) charged PAHs (0i on = 1.0) of interstellar size 
distribution (do = 3.5 A and a — 0.4; model no. 8); (b) a mixture of neutral and charged PAHs (cj>i on = 0.5, ag = 3.0 A 
and a = 0.2; model no. 9); and (c) neutral PAHs (0i on = 0, ao — 3.0 A and a = 0.2; model no. 10). Other dust 
parameters (see Table 2) are the same as those of the canonical "cold-coagulation" model (no. 1; Figure 3). 
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Fig. 9. — Size distributions for the PAH population. A log-normal functional form (see Eq.[3]) is assumed. The 
distribution is determined by a peak-parameter ao and a width-parameter a: ao = 2.5 A and a = 0.3 (solid and 
dotted lines) for the canonical best-fit model (no. 1; see Figure 3 and §3.2); ao = 3.5 A and a = 0.4 (dashed line) for 
the interstellar PAH model (see Li & Draine 2001a; see Figure 8a and §4.1 in this paper); ao = 3.0 A and a = 0.2 
(dot-dashed line) for the 0i on = 0.5 and 0i on = model (see Figures 8b,c and §4.1 in this paper). The PAH lower- 
cutoff sizes are all set at a^,? = 3.5 A except for one case (dotted line: oj^ 1 = 4.6 A [model no. 11]; see §4.2 and 
Figure 11). 
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Fig. 10. — Photodestruction time scales for PAHs at r = 20, 55, 90 AU as a function of size (upper axis) or Nc (the 
number of carbon atoms in a PAH molecule; lower axis). Shaded region depicts the age of the HD 141569A system. It 
is seen that the smallest size of surviving PAHs is ~ 4.6 A. But a stable distribution of a population of smaller PAHs 
is not impossible as long as they are continuously replenished, presumably by sublimation of icy mantles of large 
grains (in which interstellar PAHs have condensed during the dense molecular cloud phase) produced by collisions of 
large bodies such as planetesimals and comets. 
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Fig. 11. — Same as Figure 3 (model no. 1) except that the PAH component has a larger lower-cutoff size aj^ 1 = 4.6 A 
(model no. 11). In comparison with the canonical "cold-coagulation" model with = 3.5 A (model no. 1; see Figure 
3), the a^^ 1 = 4.6 A model predicts weaker 3.3 and 6.2 /xm features since this model is relatively lacking in small 
(and hence hot) PAHs (see Figure 9). 
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Fig. 12. — Upper panel (a): ratio of the radiative repulsion to the gravitational attraction (/3rp) for the best-fit 
canonical "cold-coagulation" dust (model no. 1; without ice P = 0.90 [solid line] or with ice P' = 0.73 [dashed line]). 
The long-dashed horizontal lines plot /3rp = 1,0.5 (the smaller value accounting for the fact that particles released 
at periastron from parent bodies moving on circular orbits only need to have /3rp = 0.5 for ejection [Burns et al. 
1979], assuming that there is no significant gas drag). Lower panel (b): the orbit decay timescales tpr due to the 
Poynting-Robertson drag for the best-fit canonical "cold-coagulation" dust (model no. 1; without ice P = 0.90 [solid 
line] or with ice P' — 0.73 [dashed line]) at a radial distance of r p = 203 AU from the central star (note tpr oc r 2 ; see 
Eq.[24]). The shaded region depicts the HD 141569A age (« 5 ± 3Myr). 



